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Abstract. We prove that the number of hmit cycles generated from 
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conservative perturbation of a Hamiltonian polynomial vector field 
on the plane, is bounded by a double exponential of the degree of 
the fields. This solves the long-standing infinitesimal Hilbert 16th 
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The proof uses only the fact that Abelian integrals of a given de- 
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defined over Q (the Gauss-Manin connection) with a quasiunipotent 
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1. Infinitesimal Hilbert 16th problem 

The central result of this paper is an explicit upper bound for the 
number of limit cycles born from nonsingular (smooth compact) en- 
ergy level ovals in a non-conservative polynomial perturbation of a 
polynomial Hamiltonian vector field on the plane. This problem was 
repeatedly posed in various sources under different names as the 
weakened, infinitesimal or tangential Hilbert 16th problem. In this 
introductory section we briefly outline some connections between dif- 
ferent problems concerning limit cycles of polynomial vector fields 
on the plane. Much more complete expositions can be found in the 
recent surveys [Ily02, Yak05] and the books [Zol06,IY08]. 
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1.1. Limit cycles bom by perturbations of integrable planar vector 
fields 

Limit cycles, isolated periodic (compact, nontrivial) trajectories of 
polynomial vector fields, are one of the most elusive objects of analy- 
sis. There are only a handful of tools to establish the (non)existence 
of such cycles in certain domains, all of them applying to very spe- 
cific differential equations. D. Hilbert, who included the problem of 
counting the possible number of limit cycles as the 16th problem in 
his famous list [HilOO], conjectured implicitly that the problem could 
be approached by perturbation techniques, first studying vector fields 
close to those which are "simple" from the point of view of count- 
ing their limit cycles. One such natural class is the integrable fields 
which exhibit continuous families of (nonisolated) periodic trajecto- 
ries (here and below integrability means existence of a local or global 
first integral). In [PL57] Petrovskii and Landis attempted to realize 
this program by a complexification of the problem, but their attempt 
was not successful [Ily06]. Nevertheless, the problem of estimating the 
number of limit cycles of ncar-intcgrable systems became a natural 
intermediate step towards a possible future solution of the Hilbert 
problem which still seems to be completely out of reach, see [Arn04, 
Problem 1978-6, pp. 352-363] and, most recently, [Ily08, §3.2]. 

Among the integrable systems the simplest (in many respects) 
class is that of Hamiltonian polynomial systems, vector fields cor- 
responding to a system of autonomous differential equations of the 
form 

dx dH dy dH 

di=5^^"'^^' di = -a^^"'^^' 

with a real bivariate polynomial H called the Hamiltonian; in the 
Pfaffian form these differential equations can be written as dH = 0, 
and all real level ovals of H (compact connected components of the 
level curves of the form {H = t}) are integral trajectories. 

A polynomial perturbation of the Hamiltonian system (1) can also 
be written in the Pfaffian form with a small parameter e G (M^ , 0) as 
follows, 

dH + eLO = 0, uj = P{x, y) dx + Q{x, y) dy. (2) 

In general, such perturbations destroy integrability, so that for e 7^ 
most integral trajectories will become spirals. 

We say that a (smooth) closed oval 5 C \^H = t} generates a 
limit cycle in the perturbation (2), if for any sufficiently small annu- 
lar neighborhood U oi 5 one can find arbitrarily small values of the 
parameter e such that the corresponding Pfaffian equation exhibits a 
limit cycle 6^ entirely belonging to U . If there exists a natural number 
A; ^ 1 such that for an arbitrarily narrow U and an arbitrarily small 
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£ there may coexist k limit cycles, we say that the oval 5 generates 
^ k limit cycles in the family. 

Respectively, we say that an oval 5 generates no more than k 
limit cycles in perturbation (2), k ^ 0, if there exists a small annular 
neighborhood U of 5 on the {x, y)-plane, and a small neighborhood 
V = (M^,0) of the origin on the parameter axis, such that for any 
e € y the foliation defined by the Pfaffian form dH+eto, has no more 
than k limit cycles in U . The minimal number k with this property, 
denoted by k = k{5;uj), always exists: one can easily see that for 
almost all ovals k{6;uj) = (the oval is destroyed without generating 
any limit cycle). In other words, in the sum taken over all smooth 
ovals of the real level curves {H = const} C M^, 

y{{H,u)= k{d;u)^+oo (3) 

<5C{i?=const} 

all but countably many terms mush vanish. 

It is well known that k{d; u) > only if the Poincare integral 

I = I{5,u) = j)^Pdx + Qdy, (4) 

vanishes (Poincare- Andronov-Pontryagin criterion [IY08, §26 A]^). In 
physical terms, the integral (4) is the principal asymptotic term for 
the dissipation of the energy along one period. The perturbation is 
called non-conservative if /(•, w) ^ 0. In this case a slight refinement 
of the Poincare-Andronov-Pontryagin criterion asserts that k{5;uj) 
does not exceed the multiplicity of the root of the integral I{5, u)) as 
the function of the first argument. 

The infinitesimal Hilbert problem requires to place an upper hound 
for the number of limit cycles born from nonsingular energy level 
ovals in non- conservative perturbations. The answer should depend 
only on n, i.e., the bound must be uniform over all Hamiltonians 
H of degree^ degH ^ n + 1 and polynomial 1-forms to of degree 
deg cj = max(deg P, deg Q) ^ n. Our main result solves this problem 
and gives an explicit double exponential upper bound. 

Theorem 1. The total number of limit cycles !N(//, oj) that can be 
born from all nonsingular energy level ovals of a Hamiltonian poly- 
nomial foliation in a non- conservative perturbation (2) of degree ^ n, 
is no greater than 2^^°'^^"^ . 

Here the expression Poly(n) = 0{nP) stands for an explicit poly- 
nomially growing term with the exponent p not exceeding 61. 

^ For convenience of the reader we give references to the textbook [IY08] when- 
ever possible. References to the original publications can be found in this textbook. 

^ The degrees of the Hamiltonian and the perturbation 1-form are chosen so 
that both terms in the perturbation (2) have the same degree n. 
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Besides limit cycles born from nonsingular ovals, limit cycles can 
be born from separatrix polygons (energy level curves carrying sin- 
gular points of the Hamiltonian vector field). Theorem 1 does not 
address the number of these cycles, see Remark 5 below. 



1.2. Zeros of Abelian integrals 

The Poincare integral for polynomial perturbations is an integral of 
a rational (in fact, polynomial) 1-form cu over a cycle 5 on the real 
algebraic curve {H = t}. Such integrals are called Abelian integrals 
and they can be considered as functions of all parameters occurring in 
the construction (coefficients of the 1-form and the algebraic curve). 
In particular, we can consider the Hamiltonian H and the 1-form uj 
as the parameters and look at the Poincare integral (4) as a contin- 
uous branch of a multivalued function lH,u}it) of the single variable 
t, the value taken by H on the cycle 5{t) C {H(x,y) = t}^. For a 
non-conservative perturbation, this function is not identically zero by 
definition. 

In such a context the question about limit cycles born in the per- 
turbation (2) reduces to the question about an upper bound for the 
total number of isolated zeros of real branches of the Abelian integral 
lH,Lo{t), counted with multiplicities. This upper bound should be uni- 
form over arbitrary combinations of the "parameters" H and u of 
degrees not exceeding a given natural number n. 

Theorem 1 is an immediate corollary of the following result on 
zeros of Abelian integrals. For any finite values n, m G N denote by 
?vf(n, m) the upper bound for the number of isolated real zeros of the 
integrals / counted with multiplicities, 

yi{n, m) = s\xgyi{H, uj) = sup ^ ordj Ih,u}{ • )> 

LO,H ui,H J 

degi^ ^ n -|- 1, degcj ^ m. (5) 

Here ordf 1^0 denotes the order of the root of the integral / at 
a real point t (this order is zero by definition if I{t) ^ and 1 for 
a simple root); if for a given value of t the algebraic curve {H = t} 
carries several real nonsingular ovals, the summation is extended over 
all corresponding continuous branches of I. 

Theorem 2. 

\N{n, n) ^ 2 , 

were the expression Poly(n) = O(n^) stands for an explicit polyno- 
mially growing term with the exponent p not exceeding 61. 



^ If the curve {H = t} carries several smooth real ovals, we always consider 
several branches of the integral I simultaneously. 
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This result, the first explicit uniform bound for the number of 
isolated zeros of Abelian integrals, is the most recent (hopefully, not 
the last) step in the long line of research, partly outlined in §1.4. 

Remark 1 (notation for polynomial bounds). Here and below 
we will deal with explicit bounds which involve simple or double expo- 
nentials of polynomially growing terms. In order to avoid cumbersome 
notation, we will use the following shortcuts. 

Everywhere unless explicitly waived, the symbol 0( • ) refers to 
an explicit constructive asymptotic; in particular, the notation O(n^) 
means a constructive positive function of an integer argument n which 
does not exceed the expression CnP for some explicit constant C < 
+00 and all n ^ 2. 

The notation Poly(n) stands for the constructive bound 0{nP) for 
some unspecified finite exponent p < +oo. 

Sometimes the growth rate will be estimated by the "extended 
polynomial notation" 0~^('nP), which by definition means "©(n^"*"*") 
for any positive e". A typical example is as follows: for any finite 
q>0, 

nln«n = 0+(n). (6) 

Of course, the notation 0+(nP) in the upper bound could be replaced 
by 0{nP'^^), yet when such "rounding errors" are composed, the over- 
all bound gets increased by artificial terms unrelated to the nature 
of these bounds. 

For bounds depending on several arguments, we use the notation 
0{nPm'^), Poly(n,m) and 0~^{nPm'^) in a self-explanatory way. 

The introduced notation allows to formulate the improved bounds 
in Theorems 1 and 2. What we prove in fact is the double exponential 
bound 22''°'^^"^ with Poly(n) = 0+{n^^), see (23). 

Remark 2 (accuracy of the upper bounds). The double expo- 
nential expression cited in Theorem 2 is only an upper bound with 
no claim of accuracy whatsoever. Moreover, some fine tuning of the 
tools developed in this paper, can apparently help in reducing the 
power p = 61 to a much lower value, perhaps, as low as p = 2 (at the 
price of clarity of the exposition), see §4.9. Reducing the bounds to 
less than two exponential orders would definitely require new ideas. 

On the other hand, the only known lower bounds are quadratic in 
n and linear in the degree of the form m = deg ui for a fixed n if the 
latter is allowed to grow to infinity independently of the degree of the 
Hamiltonian, cf. with (5). Thus far there is no sound conjecture on 
what might be the true asymptotic behavior of the function 3sf(n, m). 
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1.3. Bifurcations of limit cycles not covered by Theorem 1 

In order to avoid possible ambiguities, we make several remarks on what is not 
asserted in Theorem 1. The remarks below can be considered as a list of open 
problems. 

Remark 3 (Conservative and integrable perturbations). Theorem 1 gives 
no bound on the number of limit cycles if the perturbation itself is conservative, 
i.e., if the Poincare integral vanishes identically for all ovals 5 in the family. For 
instance, if lu is exact, u) = AF, F £ M.[x,y], then the entire family (2) consists of 
Hamiltonian systems and exhibits no limit cycles for all e £ (M, 0). 

On the other hand, the identical vanishing of the Poincare integral (4) in 
general does not imply that all foliations in the family (2) are integrable and 
hence do not have limit cycles at all. Indeed, the integral (4) is only the first 
variation of the Poincaxe return map with respect to the parameter e. If the 
first variation vanishes, one can compute higher variations in what is sometimes 
called the Frangoise algorithm [Fra96, Yak95], sec also [IY08, §26B]. The number 
of isolated zeros of the first not identically vanishing variation will majorize the 
number of limit cycles born from smooth ovals of the Hamiltonian field. 

For a generic polynomial H, one can show that the higher variations will 
again be Abelian integrals of certain polynomial 1-forms, yet their degrees are 
growing together with the order of the variation. Thus the problem of counting 
the limit cycles in the perturbation (2) reduces via Theorem 1 to the question on 
how many consecutive higher variations can vanish identically without the family 
(2) being necessarily integrable. This is a generalization of the famous Poincare 
center problem whose solution is unknown even in the best of the best cases 
H{x,y)=x' +y\ 

In the special degenerate (e.g., symmetric) cases the higher variations cannot 
be expressed as Abelian integrals, only by means of (polynomial expressions in) 
the so called iterated integrals [GI05,Gav05]. Formally these integrals are not 
covered by Theorem 2 below. However, L. Gavrilov and I. Iliev have shown that 
the iterated integrals still satisfy a Fuchsiam system of equations and very recently 
it was discovered that the monodromy group of this system is quasiunipotent 
[GN08]. These observations pave a way to application of Theorem 8, our principal 
result, yet many things remain to be done in order to bridge the gaps. Anyhow, 
the same question on the maximal order of the nontrivial high variations (an 
equivalent of the center problem) reappears in this case as well. 

Remark 4 (Various scenarios of integrability) . Hamiltonian vector fields 
are only the simplest case of integrable polynomial foliations on the projective 
plane RP'^. The question about all possible scenarios of integrability is one notch 
above the Poincare center problem (which addresses the question of integrability 
only locally, neax a singular point), hence is challenging and wide open. Still, some 
possibilities axe well known and documented. 

One such scenario is the Darbouxian integrability, which corresponds to re- 
placing the exact polynomial 1-form dH in (2), by a closed rational 1-form (note 
that the equations in the Pfaffian form admit multiplication by a rational factor 
without changing the behavior of the trajectories). 

Investigation of limit cycles born by perturbation of Darbouxian integrable 
systems is a completely new field, where only the first steps are now taken [Nov09, 
BM08]. One of the main difficulties is the fact that the ovals of Darbouxian inte- 
grable systems are in general not algebraic. This circumstance renders practically 
inapplicable all known tools working for Abelian integrals. In particular. Theo- 
rem 1 seems to be of no help in this context, as the Poincare integrals do not 
satisfy any finite order linear differential equation. 
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Besides Hamiltoiiiaii and Darbouxian integrability, there are some other known 
scenarios. The most difficult for analysis seems to be appearance of limit cycle by 
perturbation of symmetric systems, yet the problem is too vague to be discussed 
here. 

Remark 5 (Limit cycles born from nonsmooth level curves of the Ham- 
iltonian) . Theorem 1 asserts nothing about the number of limit cycles born from 
nonsmooth ovals, corresponding to the critical level curves of Hamiltonians. 

For a generic real Hamiltonian, the only critical level ovals arc scparatrix 
loops (liomoclinic trajectories of a nondcgcnorato saddle) and double loops (eight- 
shaped curves, butterflies), pairs of homoclinic orbits of the same saddle, which 
may generate limit cycles converging uniformly to the union of the two trajectories. 

The case of separatrix loops is well understood after the seminal works of 
R. Roussarie [Rou86,Rou89]. An upper bound for the number of cycles born from 
a simple loop can be obtained from Theorem 2, yet much weaker results on the 
maximal multiplicity of zero of Abelian integrals, in the spirit of [Mar91,Mou03] 
are sufficient and give much better bounds. A similar bound for the double loops 
could perhaps be derived using the tools from [JM94]. 

For singular level curves carrying more than one saddle critical point of H, 
one cannot in general predict the number of limit cycles based only on the first 
asymptotic terms of the Abelian integral (a substitute for the order of zero for 
points of non-analyticity of I). In [DR06,CDR07] it is shown that already for 
perturbations of a Hamiltonian foliation with two nondegenerate critical points 
on the same level curve, one can obtain limit cycles not related to zeros of Abelian 
integrals (called alien cycles). 

The general question about limit cycles born from a critical level curve of an 
arbitrary polynomial Hamiltonian, is quite challenging and essentially open. 



1.4- A few milestones 

Probably the first to realize that investigation of limit cycles in ncar-integrable 
systems may be a path to solution of the liilbert problem, were I. G. Petrovskii 
and E. M. Landis [PL55,PL57,PL59]. Although their direct approach turned out 
to be unfeasible, these seminal papers stimulated the study of perturbations of 
Hamiltonian systems. 

The first study of Abelian integrals, focused on investigation of their roots in 
connection with the bifurcation of limit cycles, was undertaken in the disserta- 
tion of Yu. Ilyashenko (adviser E. M. Landis), see [Ily69b,Ily69a]. In this work 
Ilyashenko introduced very powerful tools from complex analysis and algebraic 
geometry and implicitly formulated the infinitesimal Hilbcrt problem in the form 
we solve it now. 

One of the first explicit bounds for the number of zeros of Abelian (ellip- 
tic) integrals was obtained by R. Bogdanov [Bog76]; soon Yu. Ilyashenko gave a 
transparent proof this result [Ily78] . 

Since then the number of papers devoted to investigation of zeros of Abelian 
integrals counts in the hundreds, and it is impossible to mention even the names 
of the principal contributors. The overwhelming majority of these papers deal 
with the low-degree cases n = 3,4, where the problem is essentially settled by 
L. Gavrilov, I. Iliev and C. Li (see part 2 of the book [CL07] and references 
therein) . 

In the general case of arbitrary degree, however, very little is known. In 1984 
A. Khovanskii and A. Varclienko proved that the rmmber of zeros of Abelian 
integrals is always finite and uniformly bounded over all Hamiltonians and 1- 
forms of a given degree [Kho84, Var84], i.e., that 3Nf(n, m) < -|-oo for all finite 
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combinations n, rn. Unfortunately, the proof is purely existential and does not 
give explicit bounds on Ji{n, m) even for small degrees n, rn. 

After this celebrated result many efforts were focused on obtaining asymptotic 
constructive bounds for the "counting function" K( • ). Very soon it became clear 
that the roles of the form and the Hamiltonian are quite different from the point of 
view of the difficulty of study. More precisely, for a fixed Hamiltonian H one may 
consider forms of arbitrary degree m growing to infinity, and study the asymptotic 
of the counting function m) as m ^ +00. The first bounds for this restricted 

setting were double exponential in the degree of the form [IY95] , yet very quickly 
they were replaced by single exponential [NY95] and finally linear [Pet97] bounds. 
The ultimate result, due to Petrov and Khovanskii, is the following estimate, 

Vn, rn e N Kin, m) = OT"'{l) + 0{m). (7) 

Here the first term 05f''^'(l) is a purely existential "constant" depending on n 
(uniformly over all Hamiltonians of degree ^ n+1) and the second term 0{m) is, as 
usual, explicit and constructive. For quite some time the proof of this result existed 
only in the oral tradition, until it was published in the book [Z0IO6, Theorem 
6.26]. The proof is based on the fact that the Abelian integrals of a 1-form of 
arbitrarily high degree m over level ovals of a Hamiltonian H of degree ^ n+1 can 
always be represented as a linear combination of integrals of 1-forms of degree not 
exceeding 2n with coefficients polynomial in t of degree 0{m/n) (see Theorem 10 
below) . Based on this observation, one can conjecture that there exists an explicit 
(constructive) bound of the form 

3sf(n,m)<2^ ^* ' + 0(m) as n, m +cx) (8) 

(Theorem 2 addresses only the bound for J^{n, n)). The proof could hopefully be 
obtained by a combination of the two tecliruquos, since our methods allow placing 
an explicit upper bound on the first (existential) term in (7). 

The most recent development in connection with the infinitesimal Hilbert 16th 
problem is an explicit upper bound for the number of zeros of Abelian integrals, 
uniform over all 1-forms of degree ^ n, finite for almost all Hamiltonians H of 
degree n + 1, yet non-uniform in H. In a series of papers [Glu05,Glu06,GI06, 
GI07] A. Glutsyuk and Yu. Ilyashenko established this type of bound, which 
grows exponentially as H approaches an exceptional set of Hamiltonians. This 
result was improved in [BYIO] where a bound growing polynomial ly near the 
same exceptional set was given. 

The only class of Hamiltonians of arbitrarily high degree for which uniform 
explicit upper bounds were previously known, is the class of hyperelliptic Hamilto- 
nians of the form H{x, y) = + Q(x), Q e M.[x\. In [NY99a] it was proved that 
the number of isolated zeros of hyperelliptic integrals can be majorized by a tower 
function (iterated exponent) of n = degQ under the technical assumption that 
all critical values of the hyperelliptic Hamiltonian are real. However, the height of 
this tower was much larger than 2 (corresponding to the double exponent). 

1.5. Hyperelliptic case 

The tools developed in this paper are sufficiently flexible to give bet- 
ter results for some more specific classes of Abelian integrals. For 
instance, if instead of a general bivariatc polynomial H{x^ y) we con- 
sider only the hyperelliptic Hamiltonians of the form 

^(x,2/) = ^2/' + x"+i + Aix"-i + --- + A„, AgM", (9) 
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then some steps of the proof can be skipped or improved, see §4.9. 
As a result, we have a better bound for zeros of the corresponding 
hyperelhptic integrals. 

Theorem 3. The number of complex isolated zeros of a hyperelhptic 
Abelian integral associated with the Hamiltonian (9), is bounded by 

oO"*" (n) 

the explicit double exponent 2 

This result, completely superseding the main result of our previous 
work [NY99a], is explained in §4.9. 
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2. Background, settings, main constructions and strategy 

We begin by describing the (standard) complexification of the Abelian 
integrals and reduce the infinitesimal Hilbert 16th problem to a ques- 
tion about zeros of solutions to an integrable Pfaffian system subject 
to a condition on its monodromy. The exposition in this section pri- 
marily settles the context and notations for the main body; the recent 
textbook [IY08, §26] contains all necessary details. 

2.1. Complete Abelian integrals depending on parameters 

Let r = Fx C F'^ he the complex projective curve defined in the 
affine chart {xi,X2) on C by the equation 

H{xi,X2) = 0, H = H{x,X)= Aax" (10) 

(here and below we use the standard multiindex notation denoting by 
a G the multiindex, \a\ = ai+a2, = x'^^x'^'^). The parameters 
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A of this equation naturally vary over the complex projective space 
P*" of dimension m = ^(n + 3)(n + 2) - 1 = 0{n'^). 

For a generic combination of the parameters A the curve /\ is 
a nonsingular (smooth) Riemann surface of genus g = ^n{n — 1), 

transversal to the infinite Hne I = P^-C^. The (first) homology 
group of Fx -I in this case has the rank £ = n?, see [IY08, Theo- 
rem 26.31]. The combination of the parameters corresponding to the 
exceptional (non-smooth or non-transversal to I) curves is a pro- 
jective algebraic subset that will be denoted by Z'^,. For an arbitrary 
point A* ^ 17* one can fix a system of cycles 5i, . . . ,Si generating the 
homology Hi{rx,,'2j) with integer coefficients and transport them 
horizontally in the sense of the Gauss-Manin connexion. The result 
is a multivalued framing of the fibers Hi{rx,'L) associated with the 
topological bundle Fx ^ \ over P'^-Z'*, ramified over Z* [IY08, 
Corollary 26.28]. 

The cohomology of a generic fiber (curve) Fx is generated by 
restrictions of polynomial 1-forms on this curve. Let uja = ' 
dx2 be monomial 1-forms which are primitives of the 2-forms 
= x"dxi A dx2 with ^ ai,a2 ^ n — 1, i.e. dw^ = The 
number of such forms is exactly equal to £ = n^, and it is known 
(see Appendix A for details and references) that the LOa generate the 
cohomology of a generic fiber Fx with A ^ Z* over C. However, for 
some exceptional fibers the forms uja become linear dependent after 
restriction on Fx: the corresponding set is a proper algebraic subva- 
riety, whose union with Z"* will be denoted by Z C P"*. 

Definition 6. The period matrix X{X) is the {i x ^)-square (multi- 
valued) analytic matrix function on P"*, 



ramified over the locus Z* and nondegenerate on P"* - Z. 

Remark 7 (From projective spaces to pencils of lines) . In the formulation 
of Theorem 2 the Abelian integral occurs as a function of a distinguished variable 
t, whereas all other coefBcicnts of the Hamiltonian H are treated as parameters. In 
the definition of the period matrix all coefBcicnts of H play the same role. However, 
this difference is only superficial, and one can consider X{\) as a parametric family 
of functions of one (complex) variable, so that their isolated zeros can be counted. 

If in the expression for H in (10) all parameters A = {Xa '■ cy 7^ (0,0)} are 
fixed except for the free term Ao,o, we obtain a (complex projective) line = P^. 
These lines corresponding to different values of A G P*"^^ form a pencil of lines, 
a family of lines in P"' passing through the point in P™ with the homogeneous 
coordinates (0, . . . , 0, 1). 




(11) 
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The space P™ is biratioiially equivalent to the product P™ ^ x P'" . Moreover, 
the equivalence can be chosen to map the lines from the above pencil to the 
lines of the form Pl = {A} x P^. In suitable afHne charts, the equivalence is 

represented by the natural identification C" = C™^"' x C"', A = {X,t). After such 
identification the period matrix (11) can be indeed considered as a multivalued 
matrix function X{X, t) with singularities, defined on P^ and depending on the 
parameters A € P"*~^ which vary over a compact parameter space. 

2.2. Integrable Pfaffian systems with singularities 

Consider a smooth (nonsingular) algebraic variety M and let f2 be 
a rational (i x ^)-matrix-valued 1-form on M with a singular locus 
17 C M which is an algebraic hypersurface (eventually, itself singular 
and reducible). 

The form f2 is integrable if df2 — f2 A = 0. The integrability 
condition is necessary and sufficient for the local existence of a holo- 
morphic nondegenerate matrix solution ^( • ) for the Pfaffian system 
of equations 

dX = n-X, n £ Mat {i, (M)) (12) 

near each nonsingular point a ^ U. The local solution admits analytic 
continuation along any path 7 avoiding the singular locus. If the path 
7 is closed, the result of continuation A^X differs from the initial 
solution X by a constant nondegenerate matrix factor My G GL(£, C), 
called the monodromy associated with this path: A^X = XM^. 

Definition 8. Let r : (C\ 0) ^ (M, U), be the germ of a holomorphic 
curve, not entirely belonging to U, and a = t(0). 

A sm,all loop around the point a is a closed path which is the image 
tISs) of any sufficiently small circle = {\s\ = e} C (C, 0). Here the 
smallncss means that the image of the punctured disk {0 < |s| ^ e} 
is disjoint with U. 

A loop freely homotopic to a small loop will also be referred to 
as the small loop. All small loops "supported" by the same holo- 
morphic curve T are freely homotopic to each other. The integrabil- 
ity assumption implies that the corresponding monodromy operators 
are conjugate to each other, and in particular they have the same 
spectrum. 

Definition 9. The integrable form f2 is called quasiunipotent at a 
point a G M, if all eigenvalues of each monodromy operator associ- 
ated with any small loop around a, are roots of unity. Clearly, any 
integrable form is quasiunipotent at a nonsingular point a ^ U, since 
the corresponding monodromy is the identity. The form is (globally) 
quasiunipotent, if it is quasiunipotent at every point of the singular 
locus of f2. 
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Remark 10. If dimM = 1, i.e., the system is one-dimensional, then 
the quasiunipotence condition means that the monodromy operators 
along any sufficiently small loop around each singular point are qua- 
siunipotent. This condition can be effectively verified if the system is 
Fuchsian (has a first order pole) by inspection of the spectrum of the 
corresponding residue matrix: all eigenvalues of this matrix should 
be rational. 

However, a loop that encircles several singularities is not small, 
thus quasiunipotence does not impose any conditions on the corre- 
sponding "large" monodromy. 

The system (12) (and the corresponding matrix 1-form J?) is called 

regular on M, if for any real analytic path 7: (M,0) (M, E) the 
solution grows no faster than polynomially near the singular locus, 

\Xi7is))\^'^ ^c\s\-P VsG(]R\0), (13) 

for some real numbers c,p < +00, eventually depending on the path 
7. Analyticity of the path 7 intends to rule out spiraling along the 
singular locus. One can show that regularity is sufficient to verify 
only along (real) line segments. 

The following result can be considered as a ''^removable singularity 
assertion" for regular quasiunipotent systems. 

Theorem 4 (Kashiwara theorem [Kas81]). A regular integrable sys- 
tem that is quasiunipotent at each point outside an algebraic subset 
of codimension ^ 2, is globally quasiunipotent. 

Remark 11 (important). Ifin the definition of the quasiunipotence 
(local and global) we replace the assumption on the spectrum to con- 
sist solely of roots of unity by a weaker assumption that all these 
eigenvalues have modulus 1, then Theorem 4 remains valid as well as 
all other assertions, including our main result (Theorem 8 below). It 
is in this stronger form that the results on general Fuchsian systems 
are formulated in [BYIO]. However, we arc not aware of any natu- 
rally arising system that satisfies this weaker assumption and is not 
quasiunipotent in the usual (stronger) sense. 

Theorem 5 (folklore). The period matrix X of Abelian integrals (11) 
satisfies an integrable, regular and quasiunipotent system of equa- 
tions of the form (12) on P*". 

Proof. This is a classical "well-known fact" whose proofs are scat- 
tered over a number of classical sources. We outline only the principal 
arguments. 

The period matrix X described in (11), is monodromic: the re- 
sult of its continuation A^X along any closed path 7 avoiding the 
ramification locus 17*, differs from the initial value X by a constant 
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monodromy matrix which describes the parallel transport of the cy- 
cles 5j along the path [IY08, Corollary 26.28]. 

This implies that the logarithmic derivative Q = dX ■ is a 
single- valued matrix function defined on the complement P"* - U. 

The growth of X near the singular locus is at most polynomial: 
indeed, the length of any cycle of the integration dj{X) is growing no 
faster then polynomially along any curve and so does the integrand. 
Hence all entries of J? grow no faster than polynomially near U. Being 
single-valued, J7 is necessarily a rational matrix 1-form. The integra- 
bility condition follows immediately from the local representation of 
logarithmic derivative. 

The properties of the monodromy group of i? were studied in great 
detail. The fact that the system is quasiunipotent was proved by 
Brieskorn [Bri70] and Clemens [Cle69] by completely different meth- 
ods; the proofs were re-exposed a number of times [AGV88,Zol06]. 

The shortest way to prove the quasiunipotence is using the Kashi- 
wara theorem. A generic point of E corresponds to a simple normal 
crossing of the curve Fx or to a quadratic tangency of this curve with 
the infinite line (which is a polar locus for u). In the first case the 
monodromy along a small loop 7 around a is described by the Picard- 
Lefschetz formulas [IY08, §261]: the corresponding monodromy ma- 
trix My has a Jordan block of size 2x2 with the eigenvalue 1. In the 
second case = 1, hence all eigenvalues of are necessarily ±1. 
Both cases are clearly quasiunipotent. By the Kashiwara theorem, we 
conclude that f2 is globally quasiunipotent. 

Definition 12. Throughout this paper the Pfaffian system (12) sat- 
isfied by the period matrix (11) will be referred to as the Picard-Fuchs 
system. 

2.3. Polynomial norms 

When dealing with polynomials (both univariate and multivariate), 
we will always use the ^^-norm. 

Definition 13. The norm of a multivariate polynomial P C[zi, . . . , 
Zn], Piz) = Cq,z", Cq G C (in the standard multiindex notation) is 
the sum of absolute values of its coefficients, ||P|| = |cq,|. Clearly, 
this norm is multiplicative, 

\\PQ\\ ^ \\P\\ ■ WQl (14) 

2.4. Algebraic objects defined over Q and their complexity 

Complexification, replacing integrals over real ovals by the complex 
analytic period matrix (11), was one of the first tools of investigation 
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of Abelian integrals [Ily78]. In this section we describe an opposite (in 
a sense) step and introduce the class of differential equations defined 
over Q, the subfield of rational numbers in the field C. 

Speaking informally, an object (a polynomial, rational fraction, 
variety, Pfaffian form, differential operator, semialgebraic set etc.) 
is defined over Q, if it can be effectively constructed from the ring 
Z[zi, . . . , Zn] of lattice polynomials. For objects defined over Q, besides 
the standard algebraic characteristics like degree and dimension, one 
can always measure its size, roughly understood as the magnitude of 
the integer numbers occurring in the explicit formulas describing the 
objects. 

Perhaps the term "size" is not very successful, since the size of 
a small reciprocal 1/n G Q with n G N would be essentially equal 
to the large number n. The term "complexity" would better suit 
our purposes, but this term is too overloaded. As a compromise, the 
reader may think of the size as (the exponent of) the bitsize of the 
explicit representation of the given objects. The formal definitions 
follow. 

Definition 14. The size S {R) of a lattice (integer) polynomial P G 
Z[2;i, . . . ,Zn] is set to be equal to its norm, S (P) = ||-P||. 
The size of a rational fraction R G Q(-2;i, . . . , z„) is 

S (R) = mm{||P|| + IIQII : R = P/Q, P,Q e Z[zi, (15) 

the minimum being taken over all possible representations of -R = 
P/Q with P,Q €Z[zi,...,Zn]. 

The size of a (polynomial or rational) 1-form on defined over 
Q, is the sum of sizes of its coefficients in the standard affine chart 
C F"". 

The size of a vector or matrix rational function (resp., 1-form) 
defined over Q, is the sum of the sizes of its components. 

A parametric family of objects is defined over Q, if it is defined 
over Q on the product space x birationally equivalent to P"* 

(cf . with Remark 7) . The size of the parametric family is defined via 
the corresponding equivalence. 

One can easily continue this series of definitions, extending it 
for any class of algebraic objects and their parametric families al- 
gebraically depending on auxiliary parameters. In the future we will 
need to define the size of ordinary differential operators, see §3. 

Remark 15. The size is associated not so much with an object, but 
rather with a specific formula representing it. For instance, the poly- 
nomial 1 + t + t"^ + ■ ■ ■ + t^~^ of size n in Z[f] can be represented 
as a rational function (t" — l)/(t — 1) of size only 4 in Q(t). Yet for 
our purposes this ambiguity will not be important, since we will deal 
only with explicit upper bounds for the size which means construction 
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of formulas (representations) not involving excessively large natural 
numbers. 

The most important feature of the size of formulas defined over Q 
is its controlled increase in any algorithmically defined construction. 
For instance, the size of a sum and product of two rational fractions 
of sizes si,S2, can be estimated as follows, 

V^l 92/ V 9192 / V9l 92/ 

(16) 

Composition of maps defined over Q is also an operation that in- 
creases the size (complexity) in a limited way. 

Example 16. The projective space is birationally equivalent 

to the product P*" x P", e.g., via the standard identification C"^^"' = 
(^m ^ (^n affine charts. Such transformations do not result in 

a substantial change of the complexity of objects, in particular, the 
above equivalence does not change the complexity of rational func- 
tions defined on the corresponding birationally equivalent varieties. 

An example of effective complexity control is the following explicit 
bound on solutions of linear systems of algebraic equations. 

Example 17. Consider a parametric family of systems of linear al- 
gebraic equations of the form 

A{\)z = b{\), AgP"^, z= (zi,...,z„), (17) 

and assume that this system is defined over Q, that is, the entries 
Oij , bj of the matrix A and the right hand side vector b arc elements 
from the field Q(A) of known degrees not exceeding d, and their size 
does not exceed s. 

The system may well be non-solvable over Q(A), but in case it is 
known to have a solution, such a solution can always be found of size 
not exceeding 0(n") s". 

Indeed, after eliminating all equations that are linearly dependent 
over the field Q(A), we can represent components of some solution as 
ratios of suitable minors by the Cramer rule. Each of these minors is 
the sum of at most n! products of n entries of A(A), all of them in 
the field Q(A), so that the degree is no greater than d" and the size 
is no greater than n!(6s)" by (16). 



2.5. "Quantization principle" 



The following general principle is the primary reason why objects 
defined over Q appear in a construction which initially has no such 
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structure: Any finite bound for objects defined over Q, is explicitly 
computable in terms of their size. We give two illustrations of this 
principle. 

Example 18 (algebraic, continuation of Example 17). A non- 
homogeneous system Az = 6 of linear algebraic equations defined over 
Q may have no solution, but if the solution exists, there exists also 
a solution of norm bounded in terms of the dimension n and size s 
(complexity) of the system, \z\ ^ n!s". 

Indeed, without loss of generality we may assume that all entries 
of the matrix A and the free terms b are integer not exceeding s" 
in the absolute value. Then, as explained in Example 17, for some 
particular solution each component can be represented as the ratio of 
some minors. The numerator does not exceed n!s" (again using the 
Laplace expansion with obvious estimates), while the denominator, 
being a nonzero integer, is no smaller than 1. 

Example 19 (geometric). Assume that C M" is a basic semi- 
algebraic set defined by finitely many polynomial equalities and in- 
equalities of the form Paix) = (resp., pa ^ 0), where pa are poly- 
nomials defined over Q of degree ^ d and size no greater than s. 

The set K may well be non-compact (e.g., a half-space), but if it 
is known to be compact, its diameter can be explicitly bounded as 
follows. 

Theorem 6 ([BPR03,BV07]). If a basic semialgebraic set 

K = f]{x: pa{x) ^0}, Pa e^x], degpa ^d, \\pa\\ ^ s, (18) 

a 

is bounded, then it belongs to the ball of radius R centered at the 
origin, with R = s'^'^'"'. The constant in 0{n) is explicit. 

The same result holds for semialgebraic sets defined by polynomials 
from if s is the upper bound for their size S {pa)- 

This claim, rather obvious for the one-dimensional case n = 1, can 
be proved for arbitrary n by the algorithmic quantifier elimination 
technique (corresponding to the projection of K to M"~^ C M"). The 
quantifier elimination process can be made constructive, which results 
in a controlled increase of the complexity in each step. 

Of course, the "quantization principle" is only a guiding line, not 
a theorem; in each instance we will have to address a specific question 
on effective bounds either directly (as in Example 17) or indirectly, 
using tools of effective real algebraic geometry (as in Example 19). 

The "quantization principle" was already implicitly used in the 
proof of the general result on meandering of trajectories of arbitrary 
polynomial vector fields [NY99b]. In application to linear systems 
this principle was introduced in [Gri01,Gri03] and later in a more 
transparent and general form in [Yak;06]. 
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2. 6. Counting zeros of multivalued matrix functions of several 
variables 

The period matrix X{X), the solution of the PfafRan system (12), is 
ramified over the singular locus E. We introduce the counting func- 
tion which generalizes the number of isolated zeros of functions of 
one variable to the class of multivalued matrix-functions of several 
complex variables. 

Let = L c P"* be an arbitrary projective line not entirely 
belonging to the singular locus E of the Pfaffian system (12) on M = 
P™. The intersection L r\ S 'm. this case consists of finitely many 
isolated points, and the restriction of the matrix function X{X) on L 
will be ramified over these points. 

Let T C L-i7 be an arbitrary triangle^ an open domain bounded 
by three circular arcs eventually degenerating into line segments or 
points. Since T is simply connected, one can unambiguously choose 
a continuous holomorphic branch of the matrix function Xif) = 

{x,,{t)Y..^^,teT. 

Consider the linear span, 

e 

^x{T) = {fe ff{T) ■f=J2 ^j'^ijit)' e C}. (19) 

a finite-dimensional subspace in the space of functions of one complex 
variable (recall that T is a triangle in a complex projective line), 
holomorphic in T. 

Replacing the matrix function X{-) by a different solution X{-)M, 
M G GL(£, C) (in particular, by another branch of analytic contin- 
uation of X), does not affect the subspace ^x{T), thus the latter 
depends only on the Pfaffian matrix 1-form il. 

We define the counting function as the supremum 

7{{n)= sup sup #{t G T: /(t) = 0} ^ +CX), (20) 

taken over all triangles T disjoint with U and all nonzero functions 
from^x(r). 

Remark 20. The requirement that T is a triangle is aimed at exclud- 
ing simply connecting domains spiraling around the singular locus. 
Easy examples show that spiraling domains may contain arbitrarily 
large number of isolated zeros of very simple functions. 

Remark 21. Knowing the bounding function 3Sf(-), one can use tri- 
angulation to estimate the number of isolated roots Ji{Q/U) of any 
linear combination in any polygonal domain U C L-Z" in any line 
L E.Ito. particular, the number of real zeros ?\f(/2/ReL) can be at 
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most d ■ 'N{n), where d is the degree of f2. Indeed, the intersection 
L n U consists of no more than d points which may subdivide the 
real (projective) Hne into to no more than d intervals. Each interval 
lies inside a triangle T free from points of U, hence the number of 
isolated zeros on it does not exceed !N(J7). 

Conversely, if there is a tiling L-U = [jUihj simply connected 
polygonal domains Ui and for each domain the maximal number of 
zeros J^{Q/Ui) = ^'^'^Vfe^xiUi) #{^- /(^) — 0} finite, then one can 
easily produce an upper bound for 3sr(!7) by simply adding these 
bounds, yi{fi) ^ Y.i^{^/Ui). 

Remark 22 (semicontinuity) . When counting zeros, one can eas- 
ily pass from open to closed polygons disjoint from provided that 
the bound for the number of zeros remains uniform. 
_ Indeed, assuming that the number of zeros in any closed triangle 
T G L- U does not exceed some N, one can immediately see that the 
same bound holds also for an arbitrary open triangle T. If T contains 
A^" + 1 roots of some linear combination, one can construct a closed 
triangle T ^ T which contains all these roots, in contradiction with 
the initial assumption. 

By the same token, in the definition of the counting function (20) 
it is sufficient to consider only closed triangles T from a dense subset 

in the space of all triangles £^{F"^-E) disjoint with E. If 

sup sup #{i G T: /(i) = 0} = AT < +00, 

then 3sf(i?) is also finite and equal to N. Indeed, if some linear com- 
bination / has N + 1 roots in an "excluded" open triangle T ^ 
then one can find an arbitrarily close closed triangle T G ,5^', even- 
tually belonging to a different line L' C P™, which contains at least 
N + 1 roots counted with multiplicity in contradiction with the initial 
assumption. This follows from the fact that isolated complex roots 
of holomorphic functions cannot disappear by small perturbations by 
the Weierstrass preparation theorem. 

2. 7. Main theorem in the abstract form and discussion 

The constructive solution of the infinitesimal Hilbert 16th problem 
(Theorem 2) is obtained as a corollary to the following general fact 
about solutions of systems of differential equations. 

Theorem 7 (existential bound). Let Q he a rational matrix 1-form of 
degree d on the projective space M = P™, and (12) the corresponding 
linear system of size i x £. Assume that: 

(I) J7 is integrable; 
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(R) Q is regular; 

(U) f2 is quasiunipotent. 

Then the value of the counting function 3^{f2) is finite, 

J<{f2) < +00. 

Theorem 8 (constructive bound). In the assumptions of Theorem 7 

and the additional assumption, 

(Q) i? is defined over Q and its size is s = S 
the above finiteness is explicit: 

^ (21) 

Here Poly(d,£,m) ^ 0-^{{dtmf). 

Recall that the 0+( • )-notation was introduced in Remark 1. The 
reduction from Theorem 8 to Theorem 2 is made possible by virtue of 
the following observation which improves the "folklore" Theorem 5. 

Theorem 9 (constructive derivation). The logarithmic derivative Q = 
dX ■ of the period matrix for Abelian integrals (11), and hence 
the corresponding Picard-Fuchs system (12) is defined over Q. 

The size s = S (i?), dimension i and the degree d = degi? of the 
corresponding rational matrix function are explicitly hounded from 
above as follows, 

^^2P°iy("), d^O(n2), m^0(n2), £ = n^, (22) 

where n + 1 is the degree of the Hamiltonians. 

Proof of Theorem 2. Plugging the estimates (22) into the bound 
(21), we obtain the bound for the number of zeros of Abelian integrals, 

N(n,n) ^ 2^°Mn).2°("'^-^-^^-^) ^ 2^o(„eo,„„, ^ ^,,.,„ao,_ ^^^^ 

It remains to notice that 0+(n^°) ^ 0(n^^). This calculation illus- 
trates the need for the 0+-notation, as similar estimates will appear 
in the future. 

In fact, Theorem 7 can be relatively easily proved by the appli- 
cation of tools from the Fewnomial theory developed by A. Khovan- 
skii [Kho91] and finiteness results for analytic functions achieved by 
A. Gabrielov [Gab68]. Unfortunately, this straightforward approach 
does not allow for application of the "quantization principle" men- 
tioned in 52.5. 
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Wc give an alternative proof of Theorem 7 that admits the re- 
quired "quantization" and ultimately yields Theorem 8. The main 
ideas of this proof are outlined below in §2.8. 

The conditions of Theorem 8 are very close to optimal. Indeed, 
without the integrability assumption the system has no well-defined 
solutions. Omission of the regularity assumption allows to construct 
a linear system on (i.e., in the smallest dimension), with infinitely 
many zeros of solutions on a real interval accumulating to a singular 
point, see [Yak05]. 

The assumption of quasiunipotence also cannot be considerably 
relaxed beyond the limits indicated in Remark 11: without the as- 
sumption zeros of solutions also can accumulate to a singular point 
along the real line [Yak05]. 

2. 8. Strategy of the proof 

In this section we briefly outhne the strategy of the proof of Theorems 7, 8 and 9. 

2.8.1. The analytic core: de la Vallee Poussin theorem and its generalizations. 
The basic tool for the estimates of the number of isolated zeros is a complex analog 
of the classical de la Vallee Poussin theorem [dlVP29]. This theorem asserts that 
for a homogeneous monic linear ordinary differential equation with holomorphic 
coefficients 

j/^*' + ai(t) t/^*"'' + • ■ ■ + au{t) y^O, teycC, 

the variation of axgument of any solution y(t) along a circular arc 7 of known 
length can be explicitly bounded in terms of the uniform upper bounds Ai = 
supjg^ |ai(t)|, j = 1, . . . , fc, of the non-principal coefficients of this equation along 
the arc [Yak99, Theorem 2.6, Corollary 2.7]. This property of high order differ- 
ential equations is in stark contrast with systems of first order linear equations 
with bounded coefficients, as was discovered in [NovOl]. The assumption that the 
equation is monic (i.e., its principal coefficient is identically equal to 1) is not an 
obstruction as long as the arc does not pass near singular points of the equation, 
which correspond to zeros of the leading coefficient of a general homogeneous 
equation 

ao{t)y^''^ + --- + ak{t)y = 0. (24) 

Computing variation of argument of solutions along an arc that passes through 
(or very close to) a singular point (a root of ao( • )) is impossible in general. 

2.8.2. From a Pfaffian integrable system to an isomonodromic family of homoge- 
neous linear ordinary differential equations. The system of Pfaffian equations 
(12) can be reduced to a scalar equation of order k ^ £^ in two steps. First, we 
note that the phase space M = P*" of the system (12) is birationally equivalent 
to the product P"""^ x P^, cf. with Remark 7; moreover, this equivalence can be 
chosen so that any given projective line = L C P™ becomes a member of the 
pencil of projective lines. Denoting the parameters of the pencil by A G ^nd 
the corresponding lines by Pl = {A} x P""^, we can restrict the Pfaffian systems on 
the lines from this pencil to obtain a parametric family of Pfaffian equations on 
the projective line P^ as in Remark 7. Since the latter space is one-dimensional. 
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choosing an afRnc chart t on allows to rc- write (12) as a linear system of first or- 
der ordinary differential equations. Namely, the matrix Pfaffian 1-form restricted 
on each line pl = {A} x in the fixed chart t takes the form 

^?l{A}xpi =^Wdt, AGP™-\ (25) 
The system (12) becomes in this chart a system of linear equations 

^=A{X)-X, X = X{\,t), (26) 

with the singular loci = PI P~ (the matrix A{\, • ) is not defined if the 
entire line Pl belongs to E, yet such values of the parameter constitute a proper 
semialgebraic set in P™^^). Clearly, the regularity of the initial system (12) implies 
the regularity of all systems in the family (26). 

The condition of integrability of the initial system (either on P"' or on x 
P^, as the two are equivalent) implies that the family of the equations (26) is 
isomonodromic in the following sense. 

Let 7 C Pi^ be an arbitrary closed path in the projective line, disjoint with 
the singular locus E^^. Then by continuity there exists a small neighborhood U of 
Ao in P"""! such that for all values of the parameter from this neighborhood, the 
corresponding singular loci E^ are still disjoint with 7, and hence for all such A 
the monodromy of the system (26) along 7 is still well defined. The isomonodromy 
condition means that the corresponding operators ^ do not depend on A £ (7, 
or, more precisely, remain in the same conjugacy class. 

This condition will be especially important when the singular locus E^ un- 
dergoes a "bifurcation" at Ao, e.g., E^^ contains an isolated singularity of high 
multiplicity at to £ C (a pole of order fc ^ 2 for the corresponding matrix A(Ao, •)), 
while all nearby matrices have simple singularities (poles) at k nearby points ti(A), 
i = 1, . . . ,k. The monodromy around the circle encompassing all these points (the 
so called classical monodromy of a singularity, [AGV88]) is conjugated to the 
monodromy of a small loop around to. 

Remark 23. From now on we work only with the product space P™^^ x P"*^ and 
parametric equations and systems of equations on this product. To simplify the 
notation, we replace the parameter space P™-"^ by P*" and denote the coordinates 
on it by A instead of A. 

On the second step we reduce the parametric system of linear ordinary dif- 
ferential equations (26) to a family of scalar high order equations in the most 
straightforward way, by successive derivations and linear elimination. As a result, 
we obtain a linear equation with coefficients that are polynomials in t and the 

parameters A £ P™, 

ao(A,t):y(''' + --- + afc(A,t);y = 0, t G C, A G C". (27) 

Integrability and regularity of the initial system means that each equation in the 
family (27) is Fuchsian and the family as a whole is isomonodromic and quasi- 
unipotent. In what follows we call this family the derived equation{s). However, 
two problems arise in connection with this process: 

(1) The leading coefficient of the derived equation (27) has isolated zeros, 
in general much more numerous than the singular points of the original 
system. 



On the number of zeros of Abelian integrals 



23 



(2) For specific values of the parameters belonging to a proper algebraic subset 
5^ C P™', the leading coefficient may degenerate identically, ao(A, • ) = 
for A € ^. 

The second problem implies that as A ^ .5^ , the equation undergoes what is 
usually referred to as a singular perturbation, the situation when the coefficient 
of the highest order derivative tends to zero. Behavior of solutions of singularly 
perturbed equations may be extremely complicated, and this scenario, if it indeed 
occurs, renders the entire approach via de la Vallce Poussin's theorem unworkable, 
since after reducing to the monic form the coefficients of the corresponding linear 
equations would be large on the entire plane (or most of it), not just near singular 
points. Note that we can ignore the exceptional value of parameters \ € 5^ itself 
by virtue of the Remark 22, provided that the bound for the number of zeros 
remains uniform. 

Somewhat miraculously, the "singular perturbation" that occurs in the reduc- 
tion of a regular system to the derived equation, is only apparent: together with 
the leading coefficient, all other coefficients of the equation (27) necessarily vanish 
for the same values of the parameter ,5^ . This phenomenon was first discovered 
(in a simpler context) by A. Grigoriev [GriOl, GriOS]. 

The analysis carried out in [BYIO] is reproduced and generalized in §3 to 
show that if we consider the norms of these coefficients ||ai(A, • )|| (in the sense of 
Definition 13), then the maximum of the ratios of these norms, called the slope of 
the linear homogeneous equation, is uniformly bounded, 

max sup M^l:1|<+oo. (28) 

We note that the finiteness of this uniform bound depends crucially on the 
regularity of the original system, and does not hold in the more general context 
of [BYIO]. 

Finiteness of this supremum for a general regular family (26) implies, by the 
"quantization principle", an explicit computable bound provided that the initial 
system (26) is defined over Q. 



2.8.3. Invariant slope. The constructions described in §2.8.2, along with the 
techniques of [BYIO], provide a uniform bound for the number of zeros of so- 
lutions as long as the singular points remain well apart. A substantial difficulty 
which needs to be addressed (and cannot be circumvented for general Fuchsian 

systems) is the study of colliding singular points. 

However, under the regularity and integrability assumptions, this problem 
can be resolved. It turns out that for Fuchsian equations (operators) the slope 
as it is defined in (28) remains uniformly bounded from above even when one is 
allowed to replace the original affine chart by an arbitrary different affine (and 
more generally, even a conformally equivalent) chart on P^. This fact, somewhat 
surprising (it seems to have gone unnoticed until now), means that the collision 
of singularities can be treated by a proper resolution of singularities, as explained 
in §2.8.4. 

In Definition 29 wo introduce the notion of the invariant slope of a linear 
operator, and subsequently show that it is finite and uniformly bounded. For 
technical reasons we need to consider not only conformal changes of the indepen- 
dent variable, but also symmetrization of differential operators with respect to an 
arbitrary circle or line in P^ . The corresponding inequality is asserted in Principal 
Lemma 33, whose proof is presented in section §3. 
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2.8.4- Scale invariance and construction of an admissible configuration of slits. 
The bounds on the slope discussed in §2.8.3, imply that without loss of generality 
one may assume that the leading coefficient of the derived equation (27) is of unit 
norm with the remaining (non-leading) coefficients uniformly bounded. By the 
de la Vallee Poussin theorem, this means that the variation of argument of any 
solution can be explicitly majorized along any arc (circular or rectilinear), which 
does not pass near the singular points where the leading coefficient vanishes. 
The meaning of this dangerous proximity can bo made precise using the scale 
invariance of the invariant slope: the upper bounds for the variation of arguments 
of an arbitrary solution of the differential equation arc possible for arcs whose 
length is not very large relative to the distance to the singular locus. This scale 
invariance is one of the the key tools in the subsequent construction. For brevity 
we refer to such arcs as "short arcs" . For instance, any sufficiently small circular 
arc centered at an isolated singular point, is "short" in this sense. 

One may attempt to slit the plane with deleted singular points of the derived 
equation by such "short arcs" into finitely many simply connected domains Ui 
and apply to each domain the argument principle. This would imply an explicit 
upper bound for the counting function yi^fl/Ui), see Remark 21. Unfortunately, 
such simplistic solution is impossible, since any arc with an endpoint at a singular 
point is necessarily "long". 

To resolve this problem, we show that one can bound the number of zeros of 
(multivalued) solutions of the derived equation in punctured disks around singular 
points, and more generally in annuli, under the assumption that the monodromy 
of this equation along the equator of the annuli has eigenvalues of unit modulus 
only, and that the bounding circles are "short". Our approach goes back to the 
work [RY96] and is based on the idea called the Petrov trick after the pioneering 
work by G. Petrov [Pct90]. 

On the second step of the construction (in §4) we construct a system of "short 
arcs" subdividing the nonsingular set into simply connected domains and annuli 
bounded by nested circles. This comes in the form of a suitable clustering: we 
need to identify groups of singular points, such that distance between them is 
much smaller compared to distances to other singular points (outside the cluster). 
Then after a suitable "magnification" one cam treat points from the same clus- 
ter as "being well apart" . However, the construction needs to be iterated, since 
much smaller clusters can be parts of the larger clusters. The accurate construc- 
tion involves ideas in the spirit of the Fulton-MacPherson compactification of the 
configuration space (see [FM94]). 

The main difficulty to deal with in this second step is an effective construction 
of the system of slits so that all annuli that appear at the end will have the required 
monodromy, and the explicit calculation of the "normalized length" of these slits. 
The source of difficulty is, among other things, the apparent non-algebraicity of 
the monodromy as a function of the parameters: in general, the monodromy of 
solutions of linear systems cannot be written in closed form. The way around this 
obstacle is to use the quasiunipotence and integrability. The quasiunipotence guar- 
antees that the monodromy is quasiunipotent along the small loops which may 
encircle several colliding singularities. The integrability (isomonodromy) allows 
to conclude that the monodromy remains quasiunipotent as long as the topo- 
logical configuration of slits and singular points remains unchanged. The latter 
condition is topological (selection of a connected component in a suitable config- 
uration space), and it is known that all connected components of semialgebraic 
sets are themselves semialgebraic and effectively constructible [BPR03]. This al- 
lows the application of the "quantization principle" , transforming the existential 
finite bound for the "normalized length" of the admissible system of slits into an 
explicit upper bound for systems originally defined over Q in exactly the same 
way as was done in the first step. 
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Knowing the cxphcit length of admissible system of slits along "short" arcs 
(the "cluster diameter" as it is introduced in Definition 50) and the invariant 
slope of the equations allows to complete the proof of Theorem 8. This program 
is realized in §4. 

2.8.5. Effective derivation of the Picard-Fuchs system. To derive Theorem 2 
from Theorem 8, we need to show that the Picard-Fuchs system provided by 
Theorem 5, is in fact defined over Q for a suitably chosen parameter space. The 
arguments used in the "proof" of Theorem 5, as well as some other standard 
approaches [Gav98,Gav99], do not allow to estimate the size (complexity) of the 
Picard-Fuchs system (12). 

The necessary bounds follow from the explicit derivation of the Picard-Fuchs 
system, suggested in [Yak02] and based on an earlier work [NYOl]. We reproduce 
it below in Appendix A and derive all the required complexity estimates. 

3. From an integrable PfafRan system to an isomonodromic 
family of Fuchsian linear equations 

In this section we work with an integrable rational Pfaffian system 
(12) of dimension £ x £ and degree d on the product space P"* x P-^, 
defined over Q, of known complexity (size) s = S (i7). 

Because of the integrability, we may consider the system as an 
isomonodromic family of linear systems on P^ , parameterized by the 
parameters A G P™', and write 

dX|pi = ^2),X, Qx = -4(A, t) dt, X G P™, (29) 

in an affine chart f on P^. Denote the singular locus of the system 
(29) by C Pi ^P^ 

The main result proved in this section is the effective transforma- 
tion of the family of linear systems (29) to a parametric family of 
scalar equations of a high order, 

Dx = d'' + Ri{X,t)d''-^ + --- + RkiX,t), RjeQ{X,t), (30) 

also defined over Q, with an explicit control on the size of the family 
(defined as the sum of the sizes of all rational cocfRcicnts Rj G Q(A, t), 
cf. with Definition 14). For a family of regular operators this turns 
out to be sufficient for explicitly controlling the slope (see (28)) of 
the operators (30) uniformly over all admissible X G P"^. 

The rationale behind this step is the stark difference between sys- 
tems of first order equations and scalar high order equations in what 
concerns zeros of their solutions, see [NovOl]. To formulate the result 
more precisely, we need to elaborate the definition of the slope from 
(28) and make it conformally invariant. 
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3.1. Differential operators of higher order and their affine slope 

We will work with homogeneous linear ordinary differential equations 
with rational coefficients in the monic form 

y« + R^it) y^"-'^ + ■■■ + Rk-i{t) y' + R^{t)y = 0, 

Ri,...,Rk€C{t), (31) 

(so that the leading coefficient is identically 1) and their parametric 
families. Because of the homogeneity, the monic equation (31) can be 
re-written in the form Dy = 0, where D is a differential operator 

D = ao{t) d'' + ai{t) d''-^ + ■■■+ ak-i{t) d + ao(t), d = ^i, (32) 

at 

with polynomial coefficients ao,ai, . . . , £ C[t], ag ^ (we denote 
symbolically this fact by writing D G C [9, i] ) . Under the assumption 

gcd(ao,-- - ,ak)c[t] = 1 (33) 

the coefficients are determined uniquely modulo a scalar common 
factor. This makes the following definition self-consistent. 

Definition 24. The (affine) slope ZD of a linear ordinary differential 
operator D G C[i, 5] with polynomial coefficients as in (32), is the 
finite number 

ZD = max 1^ < +oo. (34) 
j=i,...,k \\ao\\ 

The slope of a homogeneous linear ordinary differential equation 
with rational coefficients (31) is by definition the slope of the linear 
operator D with polynomial coefficients (32)-(33), such that (31) is 
equivalent to the equation Dy = 0. 

The affine slope of an operator is a numeric measure of proximity 
of the corresponding equation to the "singular limit" , the result of 
perturbing a linear equation so that the highest derivative enters with 
a small parameter. Knowing the slope of a homogeneous equation 
suffices to place an explicit upper bound for the variation of argument 
of any nontrivial solution of this equation along an arc that does not 
pass through the singular points of this equation. 

Lemma 25 (Lemma 8 from [BYIO]). Let D be a differential operator 
(32) of order k with polynomial coefficients of degree ^ d and the slope 
S = ZD, and 7 a closed circular arc or line segment disjoint with the 
singular locus S = {ao = 0} C C, which belongs to the disk of radius 
R centered at the origin. 

Then the variation of argument of any nonzero solution of the 
homogeneous equation Dy = along the arc 7 is explicitly bounded, 

Var Arg y (t) |^ ^kS\j\ {R/rf^'^^ . (35) 
where \^\ is the length of the arc, r = dist(7, 17). 
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Remark 26. Homogeneous linear differential equations with rational 
coefficients are the natural means of describing finite dimensional lin- 
ear subspaces of holomorphic functions on P, invariant by monodromy 
around a finite locus E. For instance, if (12) is a regular Pfaffian sys- 
tem on the projective line with a singular locus S, then for any 
open set [/ C P-i7 the linear space ^x{U) spanned by components 
of any fundamental matrix solution X of the system (12) is invariant 
by the monodromy, and can be defined by the homogeneous linear 
equation 



det 



/ y xi{t) ■■■ XfXt) \ 



= 0, (36) 



where xi( •),... , Xi( • ) is a basis of ^x{U). When expanded in the 
elements of the first column and reduced to the monic form, the 
identity (36) yields a Fuchsian^ differential operator of order i with 
rational coefficients, provided that the functions Xi{t) have moderate 
growth near all points of the singular locus E [IY08, Proposition 
19.19]. 

This observation allows to define the slope of any finite-dimensional 
subspace ^ C ff{U), U C.¥-E, invariant by the monodromy trans- 
formations associated with the fundamental group 7ri(P-i7, ig); io ^ 
U , assuming that functions from ^ grow moderately near E. The 
slope is then defined as the slope of the corresponding differential 
operator D = Dj^, the differential operator with rational coefficients 
of the minimal order ord D = dime which vanishes on 

Z^ = ZD,^, Dj^€C[t,d], 
ordD^ = dimc^, V/Gif D^f = 0. ^ ' 

Note that this does not depend on the choice of operator D^, as any 
two linear differential operators with identical kernels agree up to 
multiplication of the coefficients by a common factor. The slope does 
however depend on the choice of the affine chart t. 

Remark 27. The exact choice of the simply connected domain U is 
not important as long as it is open and disjoint with the singular 
locus E, since the rational coefficients of the differential equation are 
uniquely determined by their values in any open subset of P. This 
allows us to omit the indication of the domain in the notations. 



Recall that a linear higher order differential operator is Fuchsian if it is regular, 
i.e., its solutions, multivalued functions on C Sf, exhibit polynomial growth as 
in (13). As is well-known, Fuchsian operators admit explicit description in terms 
of the order of zeros of their coefficients at the singular locus [IY08, Proposition 
19.18]. 
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3.2. Conformal invariance and symmetrization 

The notion of a slope as it is defined in (34) and (37) suffers from 
several drawbacks, the most serious being its dependence on the 
chart with respect to which the norms of the polynomial coefficients 
are computed. Applications of different conformal automorphisms 
(changes of the independent variable) of the form 

^■.t^^, detf«?Uo, (38) 



transform any linear subspace C ^{U) into another subspace 

= {ip*f: / G ^} C ^{if-'m, ^*f = fo ^. (39) 

having the same dimension and invariant by the monodromy oper- 
ators around the transformed sing ular locus S' = ip~^{S). We will 
modify the definition of the slope so that it will become invariant 
under the actions of the conformal isomorphisms (39) . 

Besides the action of conformal isomorphisms, we will need yet 
another operator on linear spaces of functions, the symmetrization (or 
conjugacy) with respect to a circle/line. The need for symmetrization 
will become clear in the context of Lemma 44. Recall that by the 
Schwarz symmetry principle, for any function / € ff{U) holomorphic 
in a domain [/ C C, the function defined in the domain by the 

formula 

f\t) = f{t), C/t = {t: t e U}, (40) 

is also holomorphic. We will refer to as the reflection of / in the 
real axis M. 

If U = is itself a domain symmetric with respect to M, then 
a C-linear subspace ^ C ff{U) will be called real (on R), if its 
reflection J^'^ = {/^ : / G ^} coincides with A finite-dimensional 
real subspace admits a basis (over C) of functions, real (i.e., taking 
real values) on i7 n M: it is sufficient to consider functions of the form 

hif + f^)- 

Any linear subspace C ^{U) can be symmetrized as follows, 

= + ^^)\unuu (41) 

(note that the functions from =Sf ® need to be restricted on the sym- 
metrization = u n u^). 

The role of the real axis can be played by any circle or real line 
7 C P equivalent to the "standard" real axis M C P by a confor- 
mal automorphism if £ Aut(P^). The conformal equivalence will be 
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denoted by the relation 7 x R. If 99 is an automorphism which trans- 
forms 7 to the real axis M, then the reflection in 7 is defined by the 

formula 

/t(<^(t)) = /(^). (42) 

Modifying all constructions above accordingly, wc arrive at the no- 
tion of a symmetrization of a linear subspace of functions Such 
symmetrization will be denoted by An subspace ^ such that 
Jf^ = =Sf is called real on 7. 

The dimension dime ■^■f' of the symmetrization depends, in gen- 
eral, on the arc 7. For instance, if a space ^ is a real on M, then 
^ = hence its symmetrization with respect to M has the same 
dimension, yet for a generic line 7^ = e'^M, < e <C 1, arbitrarily 
close to M, the dimension of the symmetrization will be twice 
the dimension of 

Definition 28. The invariant slope of a finite-dimensional monodromic 
subspace .5f C ff{U), U CF-H a. simply connected polygon, is the 
supremum of slopes of all symmetrized confornial images of L: 

<^ = supZ(^*^)®, G Aut(P^), 7xM. (43) 

Here the supremum is taken over all conformal automorphisms / and 
all circular arcs/lines 7 C P, conformally equivalent to M. 

By the natural duality between finite-dimensional monodromic 
subspaces of holomorphic functions and linear ordinary differential 
operators with rational coefficients (modulo a common factor), de- 
scribed in Remark 26, the conformal isomorphisms and symmetriza- 
tion can be defined also on differential operators with rational co- 
efficients on P^. More precisely, for a given operator D of order i 
with a singular locus 17, we denote by cp*D the operator of order 
I (unique modulo multiplication by a rational function) whose null 
space is (p*^, where =Sf = {Dy = 0} C ff{U) is the null space of D in 
any simply connected domain U C P-i7. In the same way wc denote 
by -D® the sym,m,etrization of D with respect to an arc 7 x M, the 
operator of order ^ 21 with the null space . As was already noted, 
the specific choice of the domain U is not important by Remark 27. 
An operator equal to its symmetrization, D = D®, will be called real 
on 7, though this does not mean in general that the coefficients of 
this operator (in the monic representation) indeed take real values on 
the "axis of symmetry" 7. 

Definition 29. The invariant slope <D of a linear ordinary differ- 
ential operator D with rational coefficients on is the supremum 

<D = sup Z(<^*D)©, (f G Aut(P^), 7 X M. (44) 
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Remark 30. The group Aut(P^) = PGL(2,C) of conformal automorphisms of 
the projective hne is noncompact, therefore the slope Aip* D may be unbounded 
as (p varies over this group. Similarly, the procedure of symmetrization may affect 
the slope in an uncontrollable way. 

Example 31. Let fij,: t he the linear rescaling map and Z) 6 C[3] a linear 

operator with constant coefficients. Then the slope ^^^p^D is unbounded as fj, varies 
over all nonzero complex (and even real) numbers. 

However, for Fuchsian operators the supremum in (44) is always finite. 

Proposition 32. The invariant slope of any Fuchsian operator is finite. 

We will give an indirect proof of this statement later, in Remark 39. 

The invariant slope of an operator is a semialgebraic function of 
its coefficients, yet it is very difficult to control: its value requires 
division by quantities which can be arbitrarily small. Our first main 
result circumvents this difficulty and shows that the straightforward 
reduction of a parametric linear system (29) to a parametric high 
order equation (30) results in an explicitly bounded slope. 

Principal Lemma 33. 

A. Let Q he a rational {i x i) -matrix Pfaffian system of degree d 
on P"^ X with the following properties {cf. with the assumptions of 

Theorems 7, 8), 

(I) Q is integrable; 
(R) f2 is regular. 

Then for any simply connected polygon U C F -Ex the linear 
spaces J^'x C ^{U) generated by all components Xpq{X,t) of any fun- 
damental matrix solution X of the system (12) in this domain, are 
defined by a family of linear ordinary differential equations of the 
form Dxy = 0, where D = {Dx} are Fuchsian operators as in (30). 

The family D does not depend on U and satisfies the following 
constraints: 

(i) the order k = ordD is no greater than ^ , 

(ii) the degree maxj deg is bounded by an explicit polynomial in 
£ and d = deg f2, 

(iii) the invariant slope <Dx of the operators is uniformly bounded 
over all values of the parameter A G P"* outside a proper alge- 
braic subset y ^ IP™. 

B. Under the additional assumption 

(Q) J7 is defined over Q and its size {(■om.plexity) is ,s = S (17), 
the bound for the invariant slope is explicit and double exponential, 

VA G P"^ <Dx ^ s(<^^)°^"'^ . (45) 



On the number of zeros of Abelian integrals 



31 



Note that in the proof of this result we did not use the quasiu- 
nipotence assumption. 

Definition 34. The family of equations {D\y = 0} constructed in 
Principal Lemma 33, will be referred to as the derived equation(s). 

3.3. Formal derivation 

In this section we recall a (fairly standard) reduction of the paramet- 
ric system (29) to a parametric family of linear ordinary differential 
equations, paying attention to the complexity of the algorithm. 

Lemma 35. Under the assumptions (I), (R) of the Principal Lemma 33 
all components y = Xpq{X, t) of any fundamental matrix solution X of 
the system (12) on M = x satisfy a parametric linear ordinary 
differential equation of the form D^y = as in (30) . 

The family D = {Dx}x^Frn satisfies the following constraints: 

(1) the order k = ovdD is no greater than i'^, 

(2) the degrees deg Rj do not exceed 0{i^d). 

Under the additional assumption (Q) of the Principal Lemma 33, 

(3) the family D is defined over Q and its size is hounded by a 
simple exponential, 



Proof. This claim coincides (modulo notation) with Lemma 5 from 
[BYIO], where one can find the accurate (albeit transparent) esti- 
mates. To make the exposition self-contained, we recall the main ideas 
of the proof. 

The system (12) is defined on the product space P'" x P^. Choosing an affirie 
chart t on the second factor, we can consider it as a parametric family of linear 
differential equations of the following form (cf. with (29)), 



By induction, one can instantly see that the higher order derivatives of 
the matrix X which satisfies the linear system (46), satisfy the identities 



-=M\,t)-X, Aj+i = + ■ A, j = 0,1,2,..., (47) 



where = ^i' is the identity matrix. Since A is defined over Q, all matrix functions 
in this sequence are also defined over Q. 

Derivation and matrix multiplication result in a controlled growth of the de- 
grees and sizes of the rational matrix functions Ai, A2, ■ ■ ■ : the degree grows no 
faster than linearly in j, degAj ^ jd, while the sizes S {Aj) grow no faster than 
exponentially. 



s = 



S(l?). 




(46) 
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The rational matrix {£ x i?)-functions over Q form a linear space of dimension 
f over the field Q(A, t) of rational functions on P"" x P. Thus after k f steps 
the matrices Aj will exhibit a linear dependence over this field of the form 

ylfc + RiAk-i + ■■■+ RkAo = 0, Rj € Q(A, t). (48) 

The identity (48) can be considered as a system of i'^ linear algebraic equations 
over the field Q{X,t). Solutions of this system (after elimination of all linear de- 
pendencies between equations) can be effectively computed using the Cramer rule 

as ratios of suitable determinants formed by entries of the matrices Ai, . . . ,Ak. 
This allows to place an upper bound for the sizes S (Rj) in terms of s and the 
parameters d, I. 

The differential operator D = d'' + Rk-jd\ 9 = ^ by construction van- 
ishes on X: DX = 0. Read componentwise, this matrix identity proves the Lemma. 



3.4- Existential bounds for the slope of the derived family 

The family of differential operators D = {Dx} which is constructed 

in Lemma 35, can be always reduced to a form with polynomial co- 
efficients, which will be referred to as the standard form, 

D = po(A, t) 5'^ + • • • + Pfc_i(A, t)d + pk{X, t), 
Pj G Z[A, t], gcd2[A,t](po, • • • ,Pfc) = 1, 

The operator Dx has order k for almost all values of A G P. Moreover, 
by the Bertini-Sard theorem, the subset of the parameter values A, 

^ = {X:po{X,.) = 0} 

U {A: deggcdc[t] (po(A, ■),... ,Pfe(A, • )) > 0}. (50) 

is a proper algebraic subvariety, ^ P"*, which is nowhere dense in 

pm 

For all A ^ S^, denote by a{X) = ^Dx the affine slope of the 
operator Dx- By construction, it is a semialgebraic function on the 
parameter space. 

Lemma 36. // the initial integrable system (46) is regular, the affine 
slope of the derived family Dx is globally bounded, 

sup ZDx < +00. 

pm -y 

Proof. We will prove that the function a: P"*-^ M_|_ is locally 
bounded near each point a G P"^, i.e., that there exists such neigh- 
borhood U of a, such that snpxeu-y < +00. By compactness 
of P"*, this would imply that a is globally bounded. Clearly, it is 
sufficient to consider only the points a G ,y, i.e., the values of the 
parameters A for which the leading coefficient of the operator Dx 
vanishes identically: at all other points a is continuous. 
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Moreover, we can assume without loss of generality that the pa- 
rameter A varies along a real analytic curve, 

A = 7(s), sG(M\0); 7(s)^^^s^0. 

Indeed, if the function a were not locally bounded, then the point 
(a, 0) G J7 X would be the limit point for the open semialgebraic 
set S = {(A, z) : A ^ y,a{\) > l/z}. By the curve selection lemma 
[Mil68], the point (a, 0) can be reached from inside by a real analytic 
curve s i— >• (^(s), z{s)), which means that the function a{X{s)) grows 
to infinity as s ^ 0. 

Thus wc need to consider only the particular case of a differential 
operator D G ff{R^,0)[t,d] depending on a single parameter, 

D=po{s,t)d'' + ---+pk{s,t), se{R\0), tec, d = §^, (51) 

with coefficients pj polynomial in t, real analytic in s G (]R^,0) and 
having no common factor for all s 7^ 0. Because of the real analyticity, 
we can complexify (51) to become a family of differential operators 
=const defined for all sufficiently small complex s G (C^,0). 

The singular locus of this family is the analytic curve {{s,t) G 
(C,0) X C: po{s,t) = 0}. Apart from the axis {s = 0} corresponding 
to the identically vanishing leading coefficient, each axis {s} x C inter- 
sects this curve by finitely many points corresponding to singularities 
of the equation D\^y = 0. These points lie on branches of the above 
analytic curve, therefore one can always find a disk D C C of radius 
1, such that the product (C,0) x D intersects the singular locus only 
by the disk {0} x B. 

Consider a fundamental system of solutions xi{s,t), . . . , Xk{s, t) of 
the equation Dgij = in the product (C^,0) x B. which consists of 
the linearly independent components of a matrix solution X(X(s),t) 
of the initial system (46). By the choice of B, these functions are 
holomorphic outside {s = 0} and linearly independent, but may well 
be ramified with a nontrivial monodromy A corresponding to a loop 
around the axis {s = 0}, 

A{xi, ...,Xk) = ixi,...,Xk)-M, Me GL{k, C). 

Fix a constant matrix A such that exp27riA = M (such a matrix 
always exists, since M is invertible). Then the tuple of functions 

{x[,...,x'k) = {xi,...,Xk) -s^ 

is single- valued in (C^, 0) x D. Because of the regularity of the initial 
system the new tuple of functions has at most poles of finite order 
on the axis {s = 0}. After replacing A by ^ -|- qE for sufficiently 
large g G N, we construct a tuple of functions x'j{s, t), still forming a 
fundamental system of solutions for the family {Dg} in (C, 0) x D for 
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all s 7^ 0, such that x'^ are holomorphic on the axis {s = 0}. Note that 
the restrictions of these functions on the axis itself may well become 
degenerate (linearly dependent, e.g., identical zeros). 

Consider the C-lincar subspaces Q (J^) spanned by the func- 
tions x[{ - ,s), . . . ,x'f^(- ,s). These subspaces depend holomorphically 
(in the natural sense) on s G (C^,0) as long as s 7^ 0. The above 
mentioned degeneracy theoretically means that the analyticity breaks 
down at s = 0. 

One of the keystone results, Lemma 7 from [Yak06] (cf. with 
Proposition 18.18 from [IY08]), contends that this is not the case, and 
the application s ^ is a holomorphic curve in the "Grassmanian" 
of k- dimensional subspaces in the Banach space ^(D). This result 
can be seen as a removable singularity-type theorem. One can avoid 
technical difficulties of dealing with infinite-dimensional Grassmani- 
ans by stating that one can choose different bases in the subspaces Jf^ 
which would depend analytically on s for all s G (C^,0) and remain 
linear independent as s = 0. 

Lemma 37 (Lemma 7 from [Yak06]). For any collection of holomor- 
phic functions x'i,...,x'i^ € ^((C^,0) x D) such that x'j{s, •) are 
linear independent in ^(D) for all s 7^ 0, one can construct a mero- 
morphic matrix function R{s), nondegenerate and holomorphic for 
s 7^ 0, such that the tuple of functions fj = fj{s,t), 

(/i,...,/fc) = (x'i,...,4)-i?(s) 

is holomorphic in {s, t) and linearly independent on each fiber {s = 
const}, including the axis {s = 0}. □ 

By construction, Dgfj = for all s € (C^, 0) and all j = 1, . . . ,k. 

In a standard way, one can construct a family of monic differen- 
tial operators L = {Lg} of order k with coefficients holomorphic in 
(Ci,0) xD, 

L = d''+ qi{s, t) d'''^ + ■■■+ qk-i{s, t) d + gfc(s, t), 
qjeff{{C\0)xB), j = l,...,k, 

which is also annulled by the same tuple: Lsfj = for ah s G (C\0) 
and all J = 1, . . . , fc. Since the orders of the two operators Lg and 
Dg are the same and the functions fj are linearly independent on all 
fibers {s = const}, we conclude that the operators Dg and Lg must 
be proportional, 

^4M = ^^, V.7^0, Vi = l,...,fc. (52) 
Po[s, •) 1 

Since qj are holomorphic, this implies that the ratios in the left hand 
side of (52) are holomorphic also on the axis {s = 0}. We show that 
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this implies the finiteness of the ratios of the norms in C[t\ as function 
of s ^ 0. 

Let Vj ^ be the vanishing orders of the polynomial coefficients pj 
on the axis {s = 0} : this means that pj = s'^^p'j , while pj G ^(C^ , 0) [t] 
and Pj{0, •) 0. From (52) it follows that uj ^ i/q for all j = 1, . . . , k. 
Besides, since the limit denominator p'q{0, • ) is nonzero, its norm 
\\Pq{s, •)\\c[t] is strictly positive. As a result, we conclude that the 
slope (t(s) = ZDs is continuous at s = 0: 



a{s) = max s 



Wp'jis, ■)\\c[t] 



j=i,...,k Wp'ois, -^dt]' 

This proves the local boundedness of the slope ZDg along the real 
analytic curve 7, concluding the proof of the Lemma. 



3. 5. Embedding in a conformally complete symmetric family 

The assertion of the Principal Lemma 33 concerns the invariant slope 
of the derived operator, whereas the finiteness achieved in Lemma 36 
is established only for the affine slope. Besides, bounding the invari- 
ant slope involves symmetrization, whose explicit construction may 
lead to an uncontrollable growth of the slope, sec [BYIO, Example 6]. 

To fill the gap, we embed the initial parametric family of Pfaffian 
systems (29) into a larger family of systems of larger dimension, whose 
derived equation will contain all symmetrizations of all conformal 
transforms of the initial derived family. The uniform bound for the 
affine slope in this new family gives a bound for the invariant slope 
of the original one. 

The embedding is rather straightforward in terms of the linear 
spaces. Denote by X{X,t) the matrix solution of the system (12) on 
P™ X pi and if (E Aiit(pi) ^ PGL(2, C) ^ P'^ - {a quadric} a variable 
conformal isomorphism. Then the components of the matrix function 
X{X, (f, •) = X{X, (p{-)) span for each admissible value value (A, if) G 
pm ^ p3 q£ ^j^g parameters the linear space ^Xj^ = ^*^x which is 
a conformal transform (39) of the linear space spanned by the 
components of the initial matrix function X{X, •). In other words, 
considered as a parametric family with the parameter space P™ x P^, 
the Pfaffian system with the matrix 1-form fi = dX ■ X^^ contains 
all conformal transforms of the initial system. 

The dimension and the degree of the system f2 obviously remain 
the same as that of the system f2. It is almost as easy to see that 
the size (complexity) S [fi] differs from S {f2) by involving into a 
constant power, s 1-^ s'-^^^^. Indeed, the derivatives of X with respect 
to the variables </? G P^ can be expressed over Q through entries of i?, 
t and ip using the chain rule, and it remains to apply a few times the 
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inequalities (16). Clearly, the new family is also regular, integrable 
and quasiunipotent. 

It remains to embed the family of linear spaces ^a,(^ into a larger 
family (still defined over Q albeit on a larger space) which would 
contain symmetrizations of all these spaces. Since the family ^a,<^ 
is already conformally complete (e.g., contains together with each 
space its conformal transforms), it is sufficient to symmetrize only 
with respect to a single axis, the real line M. 

The reflected matrix function X^{fi,\l^,t) = X(i2,tl>,t) is a holo- 
morphic matrix function which satisfies the rational integrable Pfaf- 
fian system dX"^ = Sl'^X'^ on x P'"^ x P^. This reflection does not 
affect neither degree, nor the dimension or the size of the Pfaffian sys- 
tem (obviously keeping it defined over Q). The block-diagonal matrix 
function 

X{\, if, n, ijj, t) = diag{X(A, t),X\n, i/;, t)} 

satisfies the integrable Pfaffian system of dimension 2i on the prod- 
uct space birationally equivalent to p2»"+6 x P^ with the coordinates 

dx = dx, ^=(^xt)' ^=(^r2t)- 

The corresponding family of subspaces ^a,i/?,^,i/) contains all sums 
■=SfA,ip+-S?J^, in particular, all symmetrizations of the conformal trans- 
forms (<^*=^a)®- 

By its explicit construction, the family (53) (considered as a Pfaf- 
fian system on p2™+6 xP^) is integrable, rational and regular. Applied 
to a system defined over Q on P"^ x P^, the construction results in a 
system again defined over Q on the larger subspace. The dimension 
is increased by the factor of 2 from i to 2i, and the size (complexity) 
is increased by an explicit constant factor 0(1). To check that the 
above doubling preserves quasiunipotence, note that a small loop (in 
the sense of Definition 8) in the product space p2™+6 x P^ projects as 
a small loop on each of the components p™+3 x P^; the correspond- 
ing monodromy is block diagonal with quasiunipotent (or identical) 
blocks. 

Thus for any Pfaffian system (12) with the matrix = {f^\} one 
can effectively construct its embedding (as a family) into a larger 
family with the Pfaffian matrix form i? = {i?^}, G p2m+6 .^^j^^jj ^j^g 
following characteristic property. 

Lemma 38. For any parameter value A € P™ there exists a param- 
eter value rj G ip2m+6 g^^/j ffiQ^f ^/jg invariant slope of the derived 
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operator Dx is equal to the affine slope of the corresponding derived 
operator Dj^, 

Thus a uniform bound for the affine slope ZDj] of the family of 

operators Dj^ is at the same time the uniform upper hound for the in- 
variant slope <Dx of the operators Dx derived from the initial family. 

3. 6. Proof of the Principal Lemma 33 

3.6.1. Proving the qualitative part. To prove the existential finite- 
ness of the latter (Part A of the Principal Lemma 33), we need to 
consider together with the initial family {f2x}, A G P™, and the cor- 
responding family Dx of derived operators the confer mal completion 
{f^T)}, V € p2"i+6^ g^mj ^j^g respective family {Dj^} as described in 
§3.5. 

Lemma 36 applied to the family {Dr]}, guarantees that the affine 
slope of these operators is bounded uniformly over 7] G p2m+6_ gy 
Lemma 38, this means that the invariant slope of the operators Dx 
is bounded by the same constant. 

3.6.2. From qualitative to quantitative bounds. It remains to prove 
Part B of the Principal Lemma and show that if the regular family 
(12) is defined over Q, then the bound for the invariant slope can be 
made explicit as follows, 

yXeFm-y <Dx ^ 8^'^^°^"'^ , s = S{D), d = degD. (54) 

Indeed, in this case the derived equation is also defined over 

Q and its size is explicitly bounded by the assertion (3) of Lemma 35 
(replacing s by its finite power s'^^^^ and i by 2i because of the differ- 
ence between the families D and D does not affect the asymptotic). 
This means that the subgraph of the affine slope function a{-) is 
a semialgebraic set defined over Q, exactly as in Example 19. Since 
the slope a{r]) = ZD^ is bounded by Lemma 36 (cf. with §3.6.1), 
Theorem 6 gives then the double exponential upper bound of the 
form s{0(i)-rf^)0(^™+'^) ^ ^(rf,)0(™) ^^^^^ a{.) on p2-+6. By 

Lemma 38, this gives the explicit uniform bound for the invariant 
slope ZDx and thus completes the proof of Principal Lemma 33. 

Remark 39 (Proof of Proposition 32). Let D be an arbitrary Fuchsian op- 
erator. Its conformal transforms (p*D, and their symmetrizations {(f*D)^ with 
respect to all arcs 7 = ip{M.) x R, (p, i/^ € Aut(P) C P^, constitute a parametric 
family Dx, A e x P^, which is a regular family of operators with compact base 
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= P*' (not necessarily defined over Q). Application of Lemma 36 to this family 
proves that the invariant slope <D is always finite; this gives a (very indirect) 
proof of Proposition 32. 

It would be interesting to achieve a direct proof together with an explicit bound 
on the invariant slope, e.g., in terms of the order of the equation, the number 
of Puchsian singularities and the absolute value of the respective characteristic 
exponents [IY08, Example 19.21]. 



4. Counting zeros of functions defined by Fuchsian 
equations 

A linear ordinary differential equation of bounded affine slope admits 

an explicit upper bound for the variation of argument of its nontrivial 
solutions along paths of bounded length, sufficiently distant from the 
singular points of the equation. For Fuchsian equations, because of 
the finiteness of the invariant slope (Proposition 32), such paths can 
be drawn with very few restrictions. One can then try and slit the 
projective line P-*^ along suitably chosen paths into "polygonal" simply 
connected domains, to which the argument principle can be applied. 
A similar approach, also based on the idea of suitable clustering, was 
suggested in [NY03] , yet its implementation there was conditioned on 
the spectral condition imposed on the monodromy group of the equa- 
tion. Unfortunately, this condition is algebraically unverifiable (one 
cannot, in general, algebraically compute the monodromy of a regu- 
lar system along an arbitrary closed loop). Moreover, this condition 
in the form required in [NY03] does not hold for the Picard-Fuchs 
system in general: there are some "large loops" whose monodromy 
does not possess the necessary spectral properties. 

In this section we suggest a way to circumvent this obstacle for 
isomonodromic families, and establish explicit upper bounds for the 
number of zeros under a weaker condition: the monodromy is required 
to be quasiunipotent only around small loops (cf. Definition 8). 



4-1. Normalized length 

We start by introducing some metric characteristics of finite config- 
urations of I points on the plane C, which are invariant under the 
action of the affine group. 

Let T be a fixed finite point subset in C. 

Definition 40. The normalized length of a closed circular arc 7 C 
C-T relative to the finite point set T C C is the finite positive 
number 
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where | • | denotes the usual Euchdean length in C and dist(-, •) the 
Euclidean distance. 

The normalized length of a line segment 7 disjoint with T is defined 
as the similar ratio 

differing only by the numeric factor 2it 6.283 .... 

The normalized length clearly is an affine invariant: for any affine 
automorphism (p G Aut(C), we have \^/T\ = \ip{'y)/ip{T)\. 

Example 41. The normalized length of all sufficiently small circles 
{\t - a\ = e}, a e T, < e ^ 1, is 1. 

The normalized length is a crude substitute for the length of a 
segment in the hyperbolic (Poincare) metric on the unit disk con- 
sidered as the universal covering space for the multiply connected 
domain C-T. However, this substitute will be more convenient for 
our purposes than the genuine hyperbolic length, among other things 
because of the semialgebraicity of the former. 

4-2. Affine invariant bounds for zeros of solutions of differential 
equations 

The notion of normalized length allows to place bounds on the vari- 
ation of argtimcnt and the mimbcr of isolated zeros for solutions of 
homogeneous ordinary differential equations in affine invariant terms. 

The following statement gives an affine invariant generalization of 
Lemma 25 in terms of the normalized length. 

Lemma 42. Let D be a Fuchsian differential operator of order k with 
coefficients of degree ^ d and the singular locus E, and 7 a closed 
circular arc or line segment disjoint with E. 

Then the variation of argument of any nonzero solution of the 
homogeneous equation Dy = along the arc 7 is explicitly bounded, 

Var Arg y{t)\^ ^ kLP^'^^ ■ <D, (56) 

where L = \^/E\ is the normalized length of the arc (see (55) and 
<D the invariant slope of the operator D. 

Proof. One can always apply an affine transformation of the affine 
plane so that the distance between 7 and the singular locus of the 
operator will be exactly 1 and the curve itself passes through the 
origin. Then the Euclidean length of the curve will not exceed L 
by the definition of the normalized length, and it will belong to the 
circle of radius L by the triangle inequality. The affine slope S of the 
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operator D in the new chart still is no greater than <D by definition 
of the invariant slope. Applying (35) with S = <D, R = \ j\ = L and 
r = 1, we obtain the inequality (56). 

In the same way as in the proof of Lemma 42, we may generalize 
the "Petrov trick" to count zeros of solutions of differential operators 
in anmili with quasiunipotcnt monodromy in affine invariant terms. 
The starting point is the following Lemma which is borrowed from 
[BYIO]. 

Lemma 43 (Lemma 10 from [BYIO]). // the monodromy of a real 
differential operator D along the equator of a symmetric annulus K = 
{p_ < \t\ < p+} has all eigenvalues on the unit circle, then the 
number of zeros of any solution in K is explicitly hounded, 

J<{D,K) ^i2k + l){2B + 1), (57) 

where k is the order of the operator and B = B{D,K) the upper 
bound for the variation of arguments of any solution of Dy = along 
the boundary circles C± = {\t\ = p±}. 

Together with Lemma 42 above, this estimate proves the following 
explicit bound for the number of zeros of solutions in annuli. 

Lemma 44. Let D be a Fuchsian operator of order k, degree ^ d and 
the singular locus S, and K C. C - E a topological annulus bounded 
by two disjoint circles C± (one of which may degenerate to a point). 

If the monodromy of D along the equator of K is quasiunipotcnt, 
then the number of isolated zeros of any solution of Dy = in K is 
explicitly bounded: 

Ji{D,K)^k^L^^'^^ ■<D, (58) 

where L is the normalized length of the boundary, L = \C-/X!\ + 
\C+/U\. 

Proof. We can always make a conformal automorphism which trans- 
forms the anmihis K into the annulus bounded by two circles cen- 
tered at the origin. Without loss of generality, replacing if neces- 
sary D by its symmetrization around the real axis R, wc may as- 
sume that D is real. The variation of argument of any solution of 
the equation Dy = along the boundary circles is bounded by 
k\C±/I!\^^'^'i<D ^ kL^^'^'i<D by Lemma 42. The inequality (57) then 
implies the bound (58). 
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4-3. Admissible systems of slits 

In this section we describe systems of arcs such that slitting the plane 
along these arcs subdivides it into components allowing for applica- 
tion of the counting tools (Lemmas 42 and 44 from §4.2). 

Definition 45. The normalized length of a union of circular arcs and 
line segments 5 = 71 U • • • U 7^, 5 C C-T, is by definition the sum 
of the normalized length of all components, 

\S/T\ = \y,/T\ + --- + \jk/T\. 

Remark 46. Note that for each term 7^ above, the normalized length 
involves the distance from the set T to 7^ and not to their union S. 
Thus the normalized length depends on the way the set S is repre- 
sented as a finite union of arcs and segments. In our constructions, 
however, this representation will always be clear from the context. 

Example 47. Let T be any two-point set and S the union of two equal circles 71,2 
centered at these points and the shortest line segment 70 connecting these circles. 

If the radii of these circles axe equal to 1/3 of the distance between the points 
of T, then the normalized length \S/T\ is equal to 3. This length can be further 
reduced to almost 2 if the radii tend to half the distance between the points. On 
the contrary, the normalized length S/T tends to infinity if the radii of the circles 
tend to 0: in this case Ijo/T] tends to infinity. 

Let, as before, T C C be a finite point set, and 5" = Sa,r C C a 
finite union of circles of the form Sa,r = ~ = '^i}) o,i ^ C, 

ri > 0, i = 1, . . . ,k. 

Definition 48. The union of circles S is called a clustering of the 
finite point set T, if all these circles are disjoint with T and pairwise 
disjoint with each other. 

A clustering subdivides points from T into nested subsets, some 
(or most) of which in principle may be empty or singletons. 
Let 5 = □■ C C - T be a clustering of T. 

Definition 49. A finite union S" = S" U 71 U • • • U 7^ of circles Si and 
line segments ji connecting them so that the circles and segments 
have only endpoints in common, is called admissible system of slits 
for a finite point set T, if the complement C-{S" U T} consists only 
of simply connected domains (of arbitrary shape) and topological 
annuli bounded by two circles (which may degenerate to a circular 
disk punctured at a point from T). 

Clearly, any clustering can be completed to an admissible system 
of slits by infinitely many ways. The number of possibilities can be 
reduced to finite, if each segment 7^ realizes the shortest path connect- 
ing the two respective circles (provided the latter are not concentric). 
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Fig. 1. Admissible system of slits around a finite point set. 

Definition 50 (principal). The cluster diameter of a finite point 

set T C C is the infimum of normahzed lengths of an admissible 
system of slits S' as in Definition 49, involving no more than a given 
number c of circular arcs: 

cdiam(r | c) = infj \S' /T\ : S' = [J Si\_\jj admissible for t|. 

By this definition, cdiam(r | c) may well be infinite, if the number 
of circular slits is too small compared to the number of points. On 
the other hand, it is obviously finite if c is sufficiently large (see the 
proof of Lemma 53) . 



4-4- Admissible system of slits for differential equations 

Let L = pQ{t)d^ + • • • + Pk-i{t)d + Pk{t) £ C[(?, t] be a differential 
operator with polynomial coefficients pi G C[t]. Its singular locus 
= {po = 0} is a finite point set which will be denoted by T. 
If S' is an admissible system of slits for the point set T = Ul in the 
sense of Definition 49, then for any topological annulus A G C-S' the 
monodromy operator M = Ma associated with the equatorial loop 
(the positively oriented loop in A which generates Tri{A)) is defined 
uniquely modulo conjugacy M i—^ C^^MC, detC ^ 0. 

Definition 51. A system of slits S' is admissible for the operator 
L, if it is admissible for the singular locus T = El and in addition 
the monodromy of L along the equatorial loop of each annulus A 
appearing in C-S" is quasiunipotent (cf. with the assumptions of 
Lemma 44). 
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The admissible system of slits for a linear system of Pfaffian equa- 
tions dx = Qxx on a projective line, is defined analogously. Obvi- 
ously, the additional constraint imposed by the requirement on the 
monodromy map is determined by the solutions, so that a system of 
slits is admissible at the same time both for a regular intcgrable ra- 
tional system i?^ on as in (29) and for the corresponding derived 
operator D\ as in (49). 

As before, we define the cluster diameter of the singular locus of 
the operator L as the infimum over all admissible systems of slits 
involving no more than c circles, 

cdiam(L | c) = infj \S'/El\ : S' = \_\Si\_\ jj admissible for l|. 

By construction, cdiam(L | c) ^ cdiam(Z'/, | c), since not all slits 
admissible for the point set T = Ul are necessarily admissible also 
for L: some anmili may have non-quasiunipotcnt monodromy. 

Let f2 be an intcgrable Pfaffian system on P™ x P and D = 
{D\}xf=pm the associated isomonodromic family of linear ordinary 
differential operators in the standard form (49), equivalent to f2 in 
the sense explained in Lemma 35. Denote, as before, by po G C[A, i] 
the leading coefficient of the family D and by y = {\: po{X, ■) = 
0} C P™" the degeneracy locus. 

Then for any X ^ ,y we have a uniquely defined finite point set 

T{\) = {t€C:po{t,X)=0}cC, X^y, (59) 

which consists of at most I = degpo ^ 0{£'^d) points (some of which 
may escape to infinity for particular values of the parameter A) . The 
linear differential operator Dx is nonsingular on C — r(A), and we can 
introduce the function 

^c(A) = cdiam(DA I c), X^y. (60) 

Principal Lemma 52. 

A. Let Q he an intcgrable rational Pfaffian system on P"* x P 
and D = {Dx} the corresponding derived isomonodromic family of 
differential operators. 

If Q is quasiunipotent and c ^ 0{i'^d), then the function #c 
troduced in (60), is globally bounded everywhere on P™. 

B. If in addition Q is defined over Q and S (1?) ^ s, then $c 
defined over Q and admits an explicit upper bound, 

^c(A) ^ VAgP--.^. (61) 

In this formulation, as before, i is the dimension of the Pfaffian system 
f2, d its degree and m the number of parameters. The proof of this 
Lemma occupies sections §4.5-§4.7. 
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4-5. Semialgebraicity of the cluster diameter 

We start by observing that the cluster diameter of an isomonodromic 
family of linear operators is (bounded by) a semialgebraic function of 
the parameters. This would be fairly easy to prove using the technique 
of quantifier elimination if the requirement on the monodromy was 
absent in the definition, since the cluster diameter of a point set is 
determined by an explicit algorithmic formula. We show that the 
isomonodromy is the key to restoring the semialgebraicity. 

Lemma 53. In the assumptions of Lemma 52A the function is 
everywhere finite on P"*-^ and semialgebraic. The subgraph of 0c 
in X is a semialgebraic set defined by polynomial equalities 
and inequalities of degree not exceeding {id)^^'^'^^ . 

In the assumptions of Lemma 52B, the function #c is defined over 
Q and has size explicitly bounded by a double exponent, i.e., its graph 
is defined by real polynomial (in) equalities with integer coefficients 
not exceeding 

S(^)^5(^<^^"'''. (62) 

Proof. Denote by I the upper bound for the degree of the coefficients 
of the derived equation D: by Lemma 35, / is bounded by 0{£'^d). 

The clusterings of the set T = {ti, ... ,ti} by c circles Si,. . . ,Sc 
can be parameterized by an open subset of the Euclidean space 

•r^^ = X = {(oi , . . . , a„ ri , . . . , re)}, (63) 

(each circle Si is defined by the equations {|i — ai| = ri}, i = 1, . . . , c). 

Consider the product P"* x The conditions ensuring that the 
circles form a clustering of the singular locus T(A) = I^Dx ^-^^ semi- 
algebraic: the singular points ti, . . . ,ti, which are algebraic functions 
of A, should satisfy the inequalities 

\ti-aj\^rj, \aj - ak\ > rj + rk or \aj - ak\ < \rj - rk\ 

for all roots ti of the leading coefficient po{t,X) and all pairs j ^ 
k, j,k = l,...,c. The latter conditions mean that the circles are 
disjoint and non-ncstcd (rcsp., disjoint and nested). Wc add to these 
conditions the inequalities ^ aj which will guarantee that the 
shortest slits connecting any two given circles, are uniquely defined 
(as no two circles arc concentric), and the conditions that the ratios 
{ti — aj)/{ti — ak) are non-real (this will guarantee that the shortest 
slits will not pass through the singular point ti). 

The points ti themselves are defined by the equation Pq{\, ti) = Q 
(roots of the leading coefficient of the differential operator Dx) and 
the degree of that leading coefficient should be maximal (equal to I) 
to avoid escape of the roots at infinity. 
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Altogether wc have 0(c + /)^ real algebraic equalities/inequalities 
of degree ^ I and size ^ s in the space of 3c + 2m = 0{c + m) real 
variables P"" x 'T^'^. 

Admissible clusterings (defined by collections of circles, without 
any reference to the monodromy yet) form a semialgebraic open sub- 
set of the total space P"* x 't^^'^, which consists of a large number of 
connected components (different topological configurations of singu- 
lar points and circles) . Each component Cq, is a semialgebraic set by 
[BPR03, Theorem 16.13]. Moreover, its algebraic complexity can be 
explicitly bounded: each connected component can be defined over 
Q by polynomial (in) equalities of degree not exceeding lO{c+m) ^^^^^ 
size (complexity) at most 5'"^*"+™' _ 

The admissible system of slits S' based on a given clustering S € 
can be achieved by a finite number of choices (the number does 
not exceed 0{c^)) of the segments jj. Indeed, one can always assume 
that the segments 7^ realize the shortest distance between each pair 
of circles of the clustering (see Fig. 1), and by construction this choice 
is uniquely defined over each connected component Ca C P"* x 
For any choice the normalized length of the resulting system will be 
a semialgebraic function on the position of the singularities and the 
clustering parameters (the coordinates of the corresponding point in 
^^'^). In other words, the normalized length of any system of slits 
|S"/r(A)| can be considered as a multivalued semialgebraic function 

#e(A, S) = {\S'/TiX)\: = 5 U 71 U • • • U 7^}, 

7j shortest slits, (64) 

on P"* X ^. Each of the finitely many continuous branches of this 
function is semialgebraic on each connected component Ca C P"^ x 'if. 

Moreover, each continuous branch of this function can be ma- 
jorized by a function defined over Q, if the initial system is de- 
fined over Q. Indeed, the distance between a point ti and the circle 
'^j — {\t ^ = fi} is Wi — \t — ai\\, the distance between two circles 
is given by a similar formula, all of them involving only the coor- 
dinate variables, the absolute value and the coefficients 0, ±1. Thus 
the normalized length of all segments is defined over Q. In the same 
manner the normalized length of the circles is defined over Q (this 
is the reason why the factor 2tt appears in Definition 40 of the nor- 
malized length). The complexity of the formula defining ^d^^S) is 
at most polynomial in m + c, since all coefficients are bounded by a 
common constant 0(1) (for brevity we denote the majorant by the 
same symbol as the initial function 4>c). 

It remains to take into account the requirement on the monodromy 
of the system (equation). In general, the monodromy of a linear equa- 
tion depends in a non-algebraic way on the equation, thus the admis- 
sibility of a system of slits cannot be defined by an algebraic condi- 
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tion. However, since the initial system is integrable, the monodromy 
is constant along any continuous branch (system of slits) based on the 
clustering varying over a connected component Ca- In other words, 
the requirement of quasiunipotence reduces to a branch selection of 
the function ^d^, S), defined on P'" x some of the branches give 
the normalized length of an admissible system of slits, while others 
do not. 

It may well happen that for a given clustering S there is no ad- 
missible system of slits based on this clustering, in which case we 
set <Pc{^,S) = +00. In the case where several systems of slits based 
on the same clustering are admissible, we can choose any of them to 
evaluate <?c(A,S'), or use the minimal value. This will not affect the 
complexity of the function ^c(A, S). 

Ultimately we can express the (major ant for the) cluster diameter 
(60) as an infimum of a semialgebraic function, 

^c(A) = inf{^c(A, S) : (A, 5) G F"^ x ^^c} ^ (65) 
s 

which itself is semialgebraic by the Taxski-Seidenberg theorem (quan- 
tifier elimination principle). Moreover, since the complexity of the 
quantifier elimination algorithm is known, we can guarantee that the 
polynomial (inequalities) defining the graph of ^c(A) over Q, have 
degree at most lO{c+m)^0{c)0{m) ^ iO{cmy' ^^^^ gj^g (complexity) at 

most . Substituting the value I = 0{i'^d), we obtain the bound 

(62). 

It remains to show that the function takes finite values for all 
values of the parameter A ^ S^, i.e., for matrices Qx (resp., operators 
-Da) with finite singular locus, provided that c^ I. 

In this case the clustering Sx which consists of exactly / circles 
centered at each singular point ti G Tx and having sufficiently small 
radius ^ 1, can be completed by finitely many segments to an ad- 
missible system of slits. Indeed, one has to connect the small disks in 
an arbitrary way with each other and with a large circle "centered at 
infinity" (the circle whose exterior contains only one singular point at 
t = 00): the only annuli that are formed by these slits, are punctured 
disks around singularities, and their monodromy is quasiunipotent 
by the quasiunipotence assumption on the initial system f? (cf. with 
Remark 10). Thus <Pc(A, Sx) < +00, hence <^c(A) < -|-oo. 

Note, however, that the finiteness of values of the function #c(A) 
for A ^ ^ docs not imply yet its local boundedness. This last step is 
achieved in the next section. 
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4.6. Local boundedness of the cluster diameter in one-parametric 
quasiunipotent families 

The cluster diameter of a finite point set T = T(A) (rcsp. a family D\ 
of equations) depending on a parameter, remains a continuous (hence 
locally bounded) function of A as long as the points of r(A) (resp. 
the singular points of -Da) do not collide. In an analytic collision 
of two or more points the topological structure of the underlying 
clusters must be chosen depending on the relative "speed" of the 
colliding points. An explicit choice of this structure is possible in 
one-parametric families. 

Lemma 54. Let Q = {i^A}Ae(ci,o) ^''^ integrable meromorphic 
quasiunipotent system on (C^,0) x and D = {D\} the correspond- 
ing derived equations with rational coefficients of degree ^ I. 

Then the function = cdiam(DA I c) is bounded over all 

A G (Ci,0), i/c > 31. 

Proof. Consider the singular locus of the system: in the coordinates 
(A, t) e (C\ 0) X C (C\ 0) X it is given by an equation p{X, t) = 
0, polynomial in t with coefficients, holomorphic on A of degree ^ /. 
Without loss of generality we may assume that p{Q, •) ^ (otherwise 
divide p by a suitable power of A) . 

The equation has / roots ti(A), . . . , t;(A), which arc algebraic func- 
tions on A and, as such, can be expressed by converging Puiseaux 
series. Passing to a fractional power of the parameter X^^^ = e, we 
may assume that each root is a holomorphic function, tj = tj(e), of 
the parameter e. The quasiunipotence of the system is preserved by 
such re-parametrization. 

We will construct a clustering Se = {(aj(e), rj(e))} G ^ of all sets 
T{e) = {ti{e), . . . ,ti{e)} for all sufficiently small values of e and an 
associated admissible system of slits for (a continuous branch of 
the function 'Pde) in the terminology of §4.5) such that the function 
<?c(£, Se) will be finite as £ — > (see (64)). 

1. The outermost circle Co of the clustering is the circle which 
contains all points of the set T(0) and is of distance at least 1 from 
them. 

The next embedded level is the union of circles Cj of radius p/2 
centered at all distinct points of the set T'(O), where 

p = mm{\ti - tj\ , tiGT(O)}. (66) 

Clearly, all these circles will be disjoint with T^e) for all sufficiently 
small £, and the normalized length relative to T{e) of these circles 
remains bounded as e ^ 0. We can add shortest slits between the 
outermost circle Co and some of the first level circles Ci to make the 
complement simply connected. 



48 



Gal Binyamini et al. 




Fig. 2. Construction of the admissible system of slits. 



Construction of the next level circles is organized in the same 
way relative to circles of the first level, so we will explain it only for 
the circle Ci around one of the points ti(0) G T(0), assuming for 
simplicity that this point is at the origin, ti{0) = 0, so that Co = 
{\t\ = l}. 

2. If among the roots tj{£) there is only one such that tj{0) = 0, 
i.e., if the origin is a "simple" (non-multiple) point of Tq, then the 
construction in this circle stops, and the (degenerate) annulus {0 < 
\t\ < 1} has finite relative length for all small e. The monodromy 
along this annulus is quasiunipotent by Remark 10. 

3. If there is more than one root tj with tj{0) = 0, then several 
holomorphic functions tj : (C^,0) —>■ C have the same 0-jet. Assume 
that these functions arc labeled as ti(e), . . . , tp(e), p ^ 2. Let k ^ 1 
be the first natural number such that fe-jets oi ti, . . . ,tp (in e) are not 
all equal between themselves. 

After the rescaling t s = {t — ti{e))/e^ in the new local chart s 
we will obtain p functions si(£:), . . . , Sp(£:), still holomorphic in £, but 
with the limits Sj(0) not all coinciding. 

Construct a circle Cg which in the chart (e, s) is large enough to 
encircle all points Sj{Q) and has distance at least 1 from them, and 
the smaller circles Cj of the form {|s — Sj(0)| = centered at each 
distinct point of the set T' = {si(0), . . . , Sp(0)}, where 

p'= min{|si-sj|, SieT'}. (67) 

In the original chart t these will be very small circles (of radius 0{e^)). 
By construction, the normalized length of Cq,C[, . . . ,C' depends 
only on the position of the points inside Ci, since all other points 
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Fig. 3. Quasiunipotence of tlie loops encircling a cluster. 



of T(e) are incomparably far. On the other hand, because the nor- 
mahzed length is affine invariant, it can be computed in the chart s, 
in which it is finite uniformly over all e ^ by the same arguments 
as in the step 1 of the proof. 

The circles Cq,C[, . . . Cp will be included in the clustering; the 
admissible system of slits is complemented by the slits between Cq 
and some of the to make the slit interior of Cq - U^=i ^'j simply 
connected (here Dj are the disks bounded by the circles Cj). As in 
step 1, these extra slits will have uniformly finite normalized length 
as e ^ 0. 

4. We need to show that the annulus bounded by Ci and Cq 
has an admissible (quasiunipotent) monodromy. This is evident if 
the fiber {e = 0} is not in the singular locus of 17, since then this 
monodromy coincides with the monodromy of Ci, which is admissible 
by the Kashiwara theorem 4. In the opposite case one needs slightly 
more involved arguments. 

Lemma 55. The annulus bounded by Ci and Cq on {e = ei}, with sufficiently 
small El, has an admissible (quasiunipotent) monodromy. 

Note that this is not an immediate corollary of the Kashiwara theorem, as the 
circle Cq has zero linking number with {e = 0}, so it cannot be boundary of a 
holomorphic disc with center at the origin. 

Proof. First, applying the translation t t — ti{e), we can assume that one of 
our curves coincides with the axis {t — 0}. 

Consider first the case where the number k which appeared on Step 3 above is 
equal to 1: this means that among the singularities forming the cluster, there are at 
least two points strictly 0(e)-distant from each other as e — > 0, \tj{e) —ti{e)\~^ = 

0(£-l). 

Consider the blow-up given in the afhne chart e 7^ by {t, e) — > (e, s = t/e), 
and denote by D = the corresponding exceptional divisor, see Fig. 3. The 
lifting of each curve t — tj{e) is the curve s = Sj{e), with curves corresponding to 
our cluster (i.e. for j — 1, . . . ,p) tending to some well-defined limit on D as e ^ 0, 
and other curves not intersecting some neighborhood U of D. The lifting of Cq is 
a circle on {e = £1} (still denoted by Cq) encompassing all points (sj(ei),ei). 
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Now, Co is homotopic along the leaf e = ci and inside t/ to a small circle 
around s = oo. Consider the second affine chart of the blow-up covering the 
neighborhood of £ = 0. In this affine chart the blow-up is given by the formulas 
{t, e) — > (t, d = s/t), the leaf {e = ei} is defined by the equation {td = ei}, and Cq 
is homotopic to the curve J = {{t = exp{2m6), d = £i exp(— 27ri^)) , 9 € [0, 1]}. 
Again, this curve has linking numbers of different signs with the t-axis and d- 
lads, so it cannot be a boundary of a holomorphic disc passing through the point 
(i,d) = (0,0), i.e. it is not a small loop. In fact, if we denote by 70 — {{t = 
exp(27ri^),d = do)} and jb = {(t = to,d = £1 oxp{2iTi9)) the two small loops 
encircling D = {t = 0} and E = {d = 0}, then 7 = jd^e^ in 7ri(F-{td = 0}), 
where F is a small neighborhood of the point (t, d) = (0, 0). 

Consider the lifting (j)* i7 of the connection O. This lifting has admissible mon- 
odromy along small loops around the strict transform of the singular locus of O. 
It also has admissible monodromy along small loops around D, since their pro- 
jections are still small loops, so their monodromy is still admissible by Kashiwara 
theorem. This in means particular that the monodromies Md and Me of 0*i7 
along jD and 7^ respectively, are quasiunipotent. But Md and Me commute 
since D and E form a normal crossing, so that tti {V \ {td = 0}) is commuta- 
tive. Thus the monodromy along 7 is equal to the product MdMe^, which is 
quasiunipotent as asserted. 

For A; > 1 one should perform k blow-ups in order to get the same situation 
near the last exceptional divisor. Again, the only curves whose strict transforms 
will intersect the last divisor D will be the curves corresponding to our cluster, 
and, deforming the loop encircling them to a neighborhood V of the point of 
intersection of D with the previous exceptional divisor, we represent it as a product 
of two monodromies. Both monodromies are quasiunipotent due to the Kashiwara 
theorem, and they commute for topological reasons, so their product is again 
quasiunipotent . 

5. One can further iterate this construction, applying it to multiple 
points of the set r'(0), if any, and constructing circles of the second 
level of embedding. Clearly, the maximal multiplicity goes down by 
at least one in each descent step, thus the tree-like clustering process 
terminates (in each branch) no later than after / steps, the total 
number of the circles in the clustering being at most 31. Thus for 
31 the conformal diameter <?c(A) is bounded over all A G (C^,0) 
as claimed. 

Remark 56. The construction of the admissible system of slits for a 
quasiunipotent integrable system is completely classical. The reader 
will easily recognize in it the "screens" of Pulton and MacPherson 
[PM94] and/or the desingularization algorithm used in the proof of 
the Kashiwara theorem [Kas81]. Definition 50 of the cluster diameter 
which is affine invariant was constructed in light of these two proofs. 

4-7. Demonstration of the Principal Lemma 52 

Consider an integrable rational family of quasiunipotent systems Q = 
{fi\\ on P"* X P-^ with the deleted indeterminacy locus ^ C P. 

The cluster diameter '^d^) is a semialgebraic function on P'^-o?^, 
continuous (hence locally bounded) outside by Lemma 53. 
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We show that this function is locally bounded at each point of 
P'" (including points of ^) for c ^ 31. Indeed, assuming that #c is 
unbounded, by the Curve Selection Lemma we can assume that 
grows to infinity along some real analytic (even algebraic) curve in 
pm_^ with an endpoint in ^ (cf. with the proof of Lemma 36). 
However, this is impossible by virtue of Lemma 54. 

Thus <P is locally bounded near each point of P"^. Because of the 
compactness of the latter, we conclude that the cluster diameter is 
globally bounded, as asserted in Part A of the Lemma. 

To prove Part B, note that by the second assertion of Lemma 53, 
the function (Pc is defined over Q by polynomial (in)equalities of de- 
gree not exceeding (£d)*^^'™^ and its complexity is explicitly bounded 

by the double exponent ^C^'^)'^*''"' as in (62), where / = 0{£'^d). By 
the "quantization principle" (Theorem 6, cf. with Example 19), the 
uniform maximum of #c does not exceed the double exponential ex- 
pression 

s I = s = s , 

which coincides with (61). □ 



4-8. Proof of Theorems 7 and 8 

Consider the Pfaffian system (12) on P"* x P^ and the corresponding 
derived family D of homogeneous differential equations (30) , and fix 
an arbitrary value of the parameters \ ^ 5^ . 

The invariant slope of the corresponding operator D\ is explic- 
itly bounded by the Principal Lemma 33, by a double exponential 
expression (45): 

<D^<,s^'^^°''^\ (68) 

In addition, the cluster diameter of the family D\ with c = 3/ = 
0{i^d) circular slits is uniformly bounded, and does not exceed the 
double exponential bound (61) 

L = cdiam(L»A I c) ^ g^o+i'^'^^)'- ^ ^gg^ 

by Principal Lemma 52. 

For each of the simply connected domains Ui formed by the admis- 
sible slits, the variation of argument does not exceed the product (56) 
(Lemma 42) which, after substitution of the bounds (68) and (69) and 

k = 0{£'^) yields the double exponential bound s^'^ ^"^^ which ab- 
sorbs all other terms in the product. This places an upper bound for 
the number of zeros 'N{f2/Ui). By the same token, the same bound 
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holds for the number of zeros J^{Q/Uj) in each annulus Uj formed by 
the admissible slits: this follows from Lemma 44. 

Since any triangle T may intersect at most 0{l) different domains 
(recall that / measures the topological complexity of the singular locus 
and hence the number of simply connected domains/ annuli, appear- 
ing by admissible slits) the total number of zeros is bounded 
by the double exponent as asserted. 

4-9. Concluding remarks 

The growth of the upper bound for L (the cluster diameter) , which is 
the fastest of the three terms in the products (56) and (58), is deter- 
mined by the complexity of the clustering algorithm. The asymptotics 
can be improved at the cost of transparency. 

For instance, among the I = 0{£'^d) singular points of the derived 
equation, almost all (except for d) are apparent singularities at which 
solutions of the equation are non-ramified (and even remain holomor- 
phic). Clearly, the location of such apparent points cannot affect the 
quasiunipotence of the monodromy along an annulus. On the other 
hand, the normalized distance from these points to the slits of an ad- 
missible system is important. Thus instead of the 31 slits used in the 
proof, one can use only 3d slits, while the corresponding normalized 
length will be given by a formula which is considerably more complex 
(involving distance between roots of an equation and a given circle). 
Yet since the complexity plays much less crucial role than the number 
of variables, the overall result will be ultimately better. 

Other more subtle modifications can be made in order to prove a 
bound which is double exponential only in the number of parameters 
of the original system. However, for the purpose of preserving the 
clarity of exposition, we shall not investigate the necessary modifica- 
tions in this paper. 

Another interesting example is that of hyperelliptic integrals. In 
this case the monodromy of any annulus free from the singular points, 
is quasiunipotent. This follows from the Lyashko-Looijenga theorem 
[Loo74] asserting that any deformation of critical values of a uni- 
variate polynomial can be achieved by a suitable deformation of its 
coefficients (a fact which is wrong for multivariate polynomials). Thus 
any k < n = degH singular points of a hyperelliptic Picard-Fuchs 
system can be isomonodromically deformed into a degenerate singu- 
larity, so any loop is homotopic to a small loop. 

This observation means that in the construction of the admissible 
system of slits one can drop the monodromy condition and hence 
skip the step of isolating a connected component Cq in the proof of 
Lemma 53. The problem of optimizing the admissible system of slits 
becomes the problem from the elementary Euclidean geometry on the 
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plane, namely the computation of the cluster diameter of a point set 
T without additional restrictions. Using elementary arguments, one 
can estimate the cluster diameter by a simple exponent of / = H^T: 
cdiam(r | 3/) ^ 2^('). 

Given that the hyperelliptic Hamiltonian of degree n + 1 has n co- 
efficients (and not O(n^), as a bivariate polynomial), we have a better 
bound for the invariant slope of the corresponding derived equation. 

oO'*' (n) 

The corresponding double exponent 2 absorbs all other depen- 

dencies and altogether the mentioned improvements give the bound 
described in Theorem 3. We leave the details to the reader. 

A. Appendix. Complexity of the Picard— Fuchs system 

In this appendix we show that the Picard-Fuchs system of linear Pfaf- 
fian equations satisfied by the periods of monomial forms is defined 
over Q and has an explicitly bounded complexity (size) and prove 
Theorem 9. 

This proof can be achieved by inspection of the effective derivation 
of the Picard-Fuchs system in [Yak02], see also [NYOl] for an earlier 
version. For the readers' convenience, we reproduce the construction 
here together with all required estimates. 

A.l. Effective decomposition in the Petrov module 

In what follows we fix a natural number n G N and denote by H the 
polynomial 



H{x,X)= Aax"GZ[x,A], x = {xi,X2), XeC^^^ 



(we use the standard multiindex notation, x = (xi, X2), a = (ai, 02) G 
For each specification of A G C"^+^ we obtain a complex poly- 
nomial i/( • , A) G C[x], which for a generic value of A is a Morse 

function with the principal homogeneous part H( - ,X) having an iso- 
lated critical point of multiplicity at the origin, where 



(the latter condition occurs if and only if i7( • , A) is square-free in 



It is well-known that the monomials x", ^ ai^2 — 1, consti- 
tute a basis in the quotient space 



0<|a|<n+l 



H{x,X)= Yl 



\a\=n+l 



C[x]). 





(70) 
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by the gradient ideal Ix for almost all (though not all) A such that 
H{ • , A) is square-free. For such values of the parameters A the forms 
jia = x'^ dxi A dx2 form the basis of the module of the Brieskorn 
lattice and their (monomial) primitives generate the so called Petrov 
module [IY08, §26E]. More precisely, we have the following result 
[IY08, Theorem 26.21]. 

Let LOa e A^{C^) be the monomial 1-forms such that dua = fJ-a, 
^ Q!i,2 ^ n — 1, and lo any other monomial form. 

Theorem 10 (see [IY08]). // the monomials x" as before generate 
the quotient space (70) for a given value X, then for any monomial 
form u there exist univariate polynomials Pa & C[t] and bivariate 
polynomials u,v E C[xi,X2] such that 

Lo = ^^(Pq o H) ■ LOa + u dH + dv ^ ai^2 ^ — 1, 

a 

dH dH 

Pa eC[t], u,v eC[xi,X2], dH = - — dxi + ^ dX2, 

oxi 0x2 (71) 

' (n + 1) degpa + degojQ 
< degv ^ degw. 

^ n + deg u 

In a similar way, the forms jia = daj„ themselves generate all 
polynomial 2'forms as a module over C[t] modulo the submodule 
dH A : any monomial 2- form fi admits a representation 

li = ^iPaoH)-fia+ dH M], r?Gyl\ (72) 

a 

with analogous inequalities between the degrees of the coefficients Pa 
and the "incomplete ratio" rj. 

While the Theorem says nothing about the dependence of the 
result of the division on A, we claim that, as functions of A, the 
polynomials Pa (the remainders) and the 1-form rj are rational and 
defined over Q if the left hand sides oj, resp., /x are over Q. For future 
calculations we will need only the situation when the expanded forms 
LO, jji are monomial of degree not exceeding 0{v?). 

Proposition 57. Assume thai lo {resp., ji) is a monomial 1-form 

{resp., 2-form) of degree at m,ost 0{n^). 

Then one can construct an expansion (71), resp., (72), so that 
the functions pa G C[t](A), the polynomials u,v € C[x, y]{X) and the 
polynomial 1-form rj € C(A) ® /l^[x,y] are all defined over Q, their 
degrees are bounded by explicit expressions growing no faster than 
0{n^) and their sizes do not exceed 2'^(" \ 
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Proof. Both systems (71), (72) arc linear with respect to the unknown 
polynomials Pa{t), u{x),v{x) and polynomial 1-form rj, thus the latter 
could all be found using the method of indeterminate coefficients. 

The degrees of these objects are explicitly bounded by Theo- 
rem 10, so we have an explicit control over the number of unknown 
indeterminate coefficients to be found, 

degpa = 0(1), dim {pq} = O(n^), dim li, dim dim ry = O(n^), 

(where by dim( • ) we mean the number of the unknown coefficients in 
the expansion of these objects as polynomials in t and x respectively). 
Altogether we see that each of the systems (71), (72) reduces to a 
system of linear (non-homogeneous) algebraic equations. The number 
of unknowns N of this system is at most 0{n'^), and all entries of the 
corresponding matrix are polynomials from Z[A]: these polynomials 
are computed using the explicit expression for H and dH respectively. 
The degrees of the entries in A are at most deg pa ^ 0(1), since H 
is linear in A by assumption. The size (complexity) of the entries is 
also bounded by n (the biggest natural number that appears in the 
expansion of dH function of A). 

Solutions of such a linear system can be obtained as in Exam- 
ple 17. By construction, the corresponding minors will be polynomi- 
als from Z[A] of degrees not exceeding 0{N) = 0{n'^) and the size 
(complexity) not exceeding N\n^ ^ 2'^("^). 

Remark 58. A more accurate analysis carried out in [Yak02] shows 
that the denominators of the rational fractions representing the poly- 
nomials Pa may involve only the parameters Aq with \a\ = n+1 cor- 
responding to the principal homogeneous part of H; the dependence 
on the non-principal coefficients with \a\ ^ n is always polynomial. 



A. 2. Gelfand-Leray derivative 

For any fixed (independent of A) polynomial 1-form lo G yl^(C^) and 
any multiindex a the derivative 

^(fuj, 6C{H = 0}, H= V A«x" (73) 

can be expressed as the integral of another (in general, only rational 
in x) 1-form 77 over the same cycle 6 C {H = 0}, if the latter satisfies 
the identity 

do; = - diJ A T]. (74) 

This follows easily [Yak02, Lemma 3] from the "standard" Gelfand- 
Leray formula which corresponds to a = (0, 0) [IY08, Theorem 26.32]. 
We use this observation to express the derivative of the period matrix 
X of the monomial 1-forms forming the basis. 
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A. 3. Effective derivation of the Picard-Fuchs system 

In this section we complete the proof of Theorem 9. 

Let X be the period i x ^-matrix of monomial 1-forms cJq, ^ 
«i,2 ^ 11 — 1, £ = n'^, cf. with (11). Denote as before lu-a = duoa- 

For any /x„ the multiple H ■ Ha \s & polynomial (in x) 2-form of 
degree ^ (n+ 1) with coefficients polynomially depending on A. 
By (72), there exist decompositions (the results of division by dH 
with remainder) 

/?-Ma= iK/3°H)-fi^+ dH Aria (75) 

0^|/3|^n+l 

in which P* = {-Pa^}a,/9 is an £ x matrix function, with entries in 
Q[t](A) (polynomial in t and rational in A). Their complexity (degree 
and size) arc bounded by Proposition 57. 

Since d{HLOa) = Hjia + dH A uia and the entries of the matrix 
function P*(t) do not depend on x, the previous identities can be 
transformed to the form 

d(Hu:a - Y,{Pa(3 oH)-UJfA =-dHA {-COa - Va), 

where d = gf^ dxi + ^ d.r2. 

Now one can choose any multiindex s € and apply the Gelfand- 
Leray formula with a = s to conclude that the partial derivative of 
the matrix function HX — P*X with respect to A^ is equal to the 
period matrix of the forms —x^{uja + rja)' 



(76) 



It remains to note the polynomial (in x) 1-forms x^{cOa + 'Va) can 
be expanded as combinations of the basic forms uja with coefficients 
in Q[i](A) of controlled degree and size by (71) and Proposition 57. 
Denote by P* = {Pai3}ap corresponding matrix functions with 
entries from Q[f](A): 

x''{uJa + Va) = 'Y{Paf3 o H) ■ UJp + Ua dH + dVa- 

Substituting these identities in (76) and integrating them over the 
cycle 5 C {H = 0}, on which the polynomial H vanishes identically, 
we conclude that 

Vs = (si, 82) G Z^, ^ si,2 ^ n - 1, ^(Po^) = Po^- (77) 
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Here Pg = P*(0) and Pq = P*(0) are matrices with entries in Q(A), 
obtained by setting t = in their initial expressions. This is the 
PfafHan integrable system on the projective space A G P"*. 

From Proposition 57 it follows that the entries of the matrix func- 
tions PqjPo are in Q(A) and have degrees growing as Poly(n) and 
the size (complexity) at most exponential (bounded by 2^°^^^'*)). This 
proves Theorem 9. 
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